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SUMMARY 


An  experimental  and  analytical  investigation  was  performed  to  characterize 
the  effects  of  underloads  on  crack  growth  and  to  evaluate  current  crack  growth  re- 
tardation models  for  predicting  underload  effects.  Experimental  data  in  the  form  of 
delay  cycles  and  overload  affected  zone  sizes  resulting  from  applications  of  a single 
overload  and  overload  combined  with  underload  were  obtained  for  approximately  200 
different  cases.  The  cases  evaluated  represented  combinations  of  load  profile  classes, 
stress  intensity  factor  parameters,  and  stress  intensity  factor  levels.  Limited  evaluations 
were  also  performed  to  determine  effects  of  hold  periods  in  tension  and  compression.  All 
experimental  work  was  performed  on  specimens  from  a single  heat  of  22I9-T85I  aluminum 
alloy  plate . 

It  was  found  that  underloads  following  overloads  decreased  delay  cycles  and  in- 
creased the  overload  shut-off  ratio.  The  decrease  was  most  pronounced  for  compressive 
underloads,  and  when  the  constant  amplitude  stress  intensity  factor  ratio,  R,  was  high. 
Underloads  proceeding  overloads  produced  essentially  the  same  results  as  the  equivalent 
tension-tension  loading  with  the  underload  truncated.  Delay  cycles  increased  as  the 
overload  increased  and  as  the  constant  amplitude  minimum  increased. 

The  overload  affected  zone  size  also  increased  as  the  overload  increased,  but 
decreased  as  the  constant  amplitude  minimum  increased.  Measured  overload  affected  zone 
sizes  were  correlated  with  overload  interaction  zone  sizes  calculated  using  the  radius  and 
diameter  assumptions  for  both  plane  strain  and  plane  stress  conditions.  Most  of  the  mea- 
sured values  were  between  the  values  computed  using  the  plane  strain  plastic  zone  radius 
and  plane  stress  plastic  zone  radius. 

Delay  cycles  increased  with  hold  time  in  tension  for  times  up  to  24  hours  which 
represented  the  longest  period  evaluated.  It  was  also  indicated  that  holding  in  tension 
may  decrease  the  overload  shut-off  ratio  for  certain  cases.  Holding  in  compression  pro- 
duced a small  reduction  in  delay  cycles. 
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Analytical  predictions  of  delay  cycles  were  made  for  all  cases  evaluated  ex- 
perimentally, Predictions  were  made  using  three  current  crack  growth  retardation 
models  in  generalized  form;  namely  the  Wheeler  model,  Willenborg  model  and  the 
Grumman  closure  model.  For  each  model,  predictions  were  made  using  an  assumed 
plastic  zone  radius  and  diameter  for  both  plane  stress  and  plane  strain  conditions. 

Based  on  correlations  of  measured  and  predicted  values,  it  was  found  that  the  best  re- 
sults were  achieved  by  pairing  models  and  overload  affected  zones  as  follows: 

1)  The  generalized  Wheeler  model  and  the  plane  strain  plastic  zone 
radius. 

2)  The  generalized  Willenborg  model  and  the  plane  stress  plastic 
zone  diameter. 

3)  The  closure  model  and  the  plane  stress  plastic  zone  radius. 

A comparison  of  the  three  models  as  paired  above  revealed  that  at  high  over- 
load ratios  both  the  generalized  Wheeler  and  closure  models  were  unconservative  while 
most  of  the  generalized  Willenborg  model  predictions  were  within  a factor  of  two  of 
measured  delay  cycles.  At  low  overload  ratios,  all  three  models  provided  reasonable 
predictions  except  for  the  tension-compression  load  class  having  high  compression  load 
for  which  the  closure  model  gave  the  best  predictions. 

Both  the  generalized  Wheeler  and  generalized  Willenborg  models  require  over- 
load shut-off  ratios  and  threshold  stress  intensity  factor  for  predictions.  Sensitivity  of  the 
models  to  these  parameters  was  evaluated.  Sensitivity  to  threshold  was  found  to  increase 
as  overload  ratio  increased,  and  sensitivity  to  shut-off  ratio  increased  as  overload  ratio 
decreased.  The  generalized  Willenborg  model  exhibited  the  least  sensitivity  to  variances 
in  both  parameters.  Predictions  were  also  made  using  three  different  crack  growth  rate 
equations  for  22I9-T85I  aluminum  alloy.  Generally,  there  was  little  difference  between 
predictions  for  the  different  equations. 
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SECTION  I 
INTRODUCTION 


Analytical  methodology  to  predict  subcritical  crack  growth  in  aerospace  struc- 
tures subjected  to  complex  loading  is  an  essential  element  in  the  overall  fracture  control 
program  currently  being  applied  on  fracture  critical  structures.  Current  analytical  methods 
for  complex  loading  are  not  precise;  however,  considerable  effort  has,  and  is,  being 
directed  toward  a better  understanding  of  the  crack  growth  process.  From  such  efforts  will 
evolve  improved  analytical  modeling  of  fatigue  crack  behavior.  Within  the  current  state 
of  technology,  an  unconservative  or  conservative  crack  growth  prediction  is  possible  de- 
pending upon  the  analysis  methodology  selected  for  a particular  load  profile.  Conse- 
quently, current  methods  must  be  judiciously  applied  and  substantiated  by  adequate  test- 
ing. Eventually,  unconservatism  must  be  eliminated  for  reasons  of  safety  and  structural 
life,  and  over-conservatism  must  be  eliminated  since  it  integrates  throughout  the  entire 
design  process  to  adversely  effect  total  performance. 

It  is  generally  agreed  that  linear  analysis  produces  acceptable  crack  growth  pre- 
dictions for  constant  amplitude  loading  provided  an  adequate  data  base  is  available,  a 
valid  stress  intensity  factor  solution  is  available,  and  environmental  variations  are  ex- 
cluded. The  introduction  of  a high  load,  however,  retards  subsequent  crack  growth  to 
the  extent  that  use  of  linear  analysis  without  considering  load  interaction  effects  is  pre- 
cluded. Considerable  data  are  available  characterizing  this  retardation  effect,  and 
several  retardation  models  have  been  developed  for  predicting  the  overload  effects  on 
subsequent  crack  growth  rate.  Limited  data  are  available  which  show  that  this  retarda- 
tion effect  is  reduced  when  the  high  load  is  followed  by  a compressive  load  or  underload. 
The  current  retardation  models  account  for  the  beneficial  effects  produced  by  overload 
conditions,  but  most  models  do  not  consider  a reduction  in  these  beneficial  effects  when 
underloads  are  included.  Further  characterization  of  the  underload  effect  and  its 
interaction  with  overloads  is  necessary  to  identify  controlling  parameters  which  will  guide 
future  development  of  analytical  methods  to  more  accurately  predict  crack  growth  under 
complex  loading. 

This  program  represents  an  experimental  and  analytical  investigation  to  charac- 
terize crack  growth  behavior  associated  with  underloads  and  their  interaction  with  over- 


loads.  The  experimental  effort  encompassed  approximately  200  variables  which  rep- 
resented combinations  of  load  profile  classes  and  stress  intensity  factor  parameters. 

All  experimental  evaluations  were  performed  on  specimens  from  a single  heat  of  2219- 
T85I  aluminum  alloy  plate.  Analytical  predictions  were  made  for  all  variables  evalua- 
ted experimentally.  Predictions  were  made  using  three  current  retardation  models  in 
generalized  form,  and  these  predictions  were  correlated  with  the  experimental  data. 
Sensitivities  of  the  models  to  such  parameters  as  overload  shut-off  ratio,  stress  intensity 
factor  threshold  and  overload  affected  zone  size  formulations  were  evaluated. 


2 


SECTION  II 
SCOPE 


1.  LOAD  PROFILES 

An  experimental  program  was  performed  to  characterize  effects  of  underloads 
on  crack  growth  and  to  provide  data  with  which  to  correlate  retardation  model  predic- 
tions. Experimental  data  in  the  form  of  delay  cycles  and  overload  affected  zone  sizes 
resulting  from  application  of  a single  overload,  underload,  and  overload  combined  with 
underload  were  obtained  for  approximately  200  different  cases.  The  cases  evaluated 
represented  combinations  of  load  profile  classes,  stress  intensity  factor  parameters,  and 
stress  intensity  factor  levels.  The  load  profile  classes  evaluated  are  defined  in  Figure  I 
and  include  (a)  tension-tension;  (b)  tension-zero;  (c)  zero-tension;  (d)  tension-compression; 
and  (e)  compression-tension.  Although  the  retardation  models  currently  do  not  account 
for  hold  time  effects,  classes  (f)  and  (g)  were  included  to  experimentally  evaluate  the 
significance  of  hold  periods.  These  load  .classes  were  illustrated  in  terms  of  stress  in- 
tensity factor,  K,  since  the  experimental  evaluations  were  performed  under  quasi-constant 
K conditions  as  described  in  a subsequent  section. 

2 . STRESS  i NTENSITY  FACTOR  PARAMETERS 

Generally,  the  program  was  designed  to  systematically  determine  the  effects  of 
the  following  parameters  on  subsequent  crack  growth  for  each  of  the  load  profile  classes: 

(a)  the  overload  stress  intensity  factor  ratio,  S-Kj/K2;  (b)  the  constant  amplitude  stress 
intensity  factor  ratio,  maximum  of  the  constant  amplitude  cycle,  K2; 

(d)  the  underload  minimum,  or  K^}  and  (e)  the  overload-underload  and  underload- 
overload  sequences.  The  program  was  also  structured  to  determine  overload  shut-off 
ratios  and  the  effect  of  underload  on  shut-off  ratio.  The  stress  intensity  factor  parameters 
and  levels  referenced  above  are  defined  in  Figure  2. 

3.  BASIC  PROGRAM 

The  basic  experimental  program  was  established  to  provide  data  for  all  parame- 
ters of  interest  except  variations  in  K2  and  hold  time  effects.  As  shown  in  Table  I all 
five  load  profile  classes  were  included  which  allowed  evaluations  of  underload  magnitude 
effects  as  well  as  the  overload-underload  and  underload-overload  sequences.  I<2  was 
held  constant  at  10  KSI\/Tn.  for  the  basic  program,  and  that  value  was  selected  to  pro- 
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vide  growth  rates  In  the  10  to  10  inch  per  cycle  range  for  the  constant  amplitude 
conditions  contained  in  the  basic  program.  Several  values  of  overload  stress  intensity 
factor  ratio,  S,  were  included  for  each  of  three  different  constant  amplitude  stress  in- 
tensity factor  ratios,  R=0.l,  0.3  and  0.5.  The  number  of  S values  for  each  combination 
of  load  profile  class  and  R depended  upon  the  overload  shut-off  ratio  as  later  discussed 
under  TESTING  PROCEDURES.  With  K2  constant,  the  different  R values  provided  varia- 
tions in  AK  and  K^.  Conditions  ranging  from  overload  without  underload  to  overload 
coupled  with  high  compressive  underload  were  provided  by  varying  the  underload  ratio, 

U,  and  compressive  underload  ratio,  U^. 

4.  SUPPLEMENTARY  PROGRAM 

In  order  to  determine  K2  effects,  the  basic  program  was  supplemented  by  the 
K2“7.78  and  14  KSiyfnT  evaluations  defined  in  Tables  II  and  III,  respectively.  A dif- 
ferent R value  was  assigned  to  each  K2  such  that  AK  was  constant  at  7 KSIv/lru  When 
combined  with  the  K2-IO  KSiyiii,  cases  for  R=0.3  in  the  basic  program,  three  different 
values  of  K2  resulted  for  AK  constant  at  7 KSI \/in . Several  values  of  overload  stress  in- 
tensity factor  ratio  were  included  as  was  done  in  the  basic  program.  Only  four  of  the  five 
load  profile  classes  were  included.  The  zero-tension  class  was  eliminated  since  results 
from  the  basic  program  were  essentially  the  same  for  tension-tension  without  underload, 
zero-tension  and  compression-tension.  The  more  severe  compression-tension  load  pro- 
file class,  however,  was  retained  to  evaluate  underload  proceeding  overload  for  the 
different  values  of  K2. 

5.  HOLD  PERIOD  PROGRAM 

As  shown  in  Table  IV,  the  tension-tension  load  profile  class  was  employed  to 
evaluate  effects  of  hold  periods  in  tension  at  Kj  for  times  up  to  24  hours.  Similarly, 
the  tension-compression  load  profile  class  was  employed  to  evaluate  effects  of  holding 
in  compression  at  K^.  The  stress  intensity  factor  parameters  for  each  load  profile  class 
were  selected  to  provide  two  different  values  for  each  stress  intensity  factor,  K,  and  to 
represent  cases  below  the  overload  shut-off  ratio  without  hold  time  effects. 
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SECTION  III 

EXPERIMENTAL  PROGRAM 

1 . MATERIAL  AND  SPECIMENS 

xAll  experimental  evaluations  were  performed  on  plates  from  a single  heat  of 
22I9-‘T851  aluminum  alloy.  The  plates  were  from  the  same  heat  as  those  used  on  a 
previous  program,  Reference  I,  to  evaluate  crack  growth  for  arbitrary  spectrum  loading. 
Basic  material  characterization  tests  were  performed  as  part  of  the  previous  program,  and 
the  results  are  documented  in  Reference  I . The  plates  were  48  inches  wide  by  72  inches 
long  and  had  a nominal  thickness  of  5/8-inch. 

Test  specimen  blanks  were  machined  from  six  plates,  and  each  blank  was  uniquely 
serialized  to  identify  the  plate  from  which  it  came  and  its  location  within  the  plate  as 
shown  in  Figure  3.  The  blanks  were  then  machined  to  produce  six  inch  wide  center  crack 
specimens  of  the  configuration  shown  in  Figure  4.  AH  test  section  thicknesses  were  0.250 
inch,  and  the  longitudinal  grain  was  oriented  parallel  to  the  loading  direction.  A 0.5 
inch  long  crack  starter  slot  was  placed  in  the  center  of  each  specimen  by  electrical- 
discharge-machining  (EDM). 

2.  TESTING  PROCEDURES 

All  testing  was  performed  in  two  identical  electrohydraulic  servo  controlled 
test  systems.  Each  system  contained  the  necessary  electronic  elements  properly  integra- 
ted to  provide  control  of  the  servo  loop,  monitor  loads  and  perform  failsafe  functions. 
Additionally,  each  system  was  interfaced  to  a digital  computer  which  was  used  to  program 
the  tests  as  well  as  store  and  reduce  data.  For  tests  requiring  compressive  loads,  the 
specimens  were  fitted  with  Teflon  lined  support  bars  to  prevent  buckling. 

Crack  length  measurements  were  made  on  one  surface  of  the  specimens  using  an 
optical  technique.  Measurements  were  made  to  within  ± 0.001  inch  by  microscopically 
observing  the  crack  tip  against  a reference  grill  lightly  engraved  on  the  specimen  surface. 
Measurements  were  made  while  holding  at  the  constant  amplitude  mean  for  the  particular 
test  being  performed. 
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Load  form  was  sine  wave  in  all  cases,  and  all  constant  amplitude  cycles  were 
applied  at  a rate  of  12  cycles  per  second.  Overloads  were  applied  at  one-quarter  cycle 
per  second  as  well  as  were  underloads.  The  test  area  of  each  specimen  was  enclosed  in 
a clear  acrylic  chamber  containing  dessicant,  and  testing  was  subsequently  performed  at 
ambient  temperature,  70-80  degrees  F . 

The  EDM  crack  starter  slot  in  each  specimen  was  precracked  to  produce  an  ini- 
tial crack  length,  2a,  of  1.0  inch.  Loads  were  continually  reduced  during  precracking, 
and  the  final  loads  represented  the  constant  amplitude  condition  for  the  particular  test 
to  be  performed . 

3.  TEST  SYSTEM  COMPUTER  PROGRAM 

Each  test  system  computer  was  programmed  to  accept  desired  stress  intensity 
factor  values  as  input,  compute  corresponding  loads  using  the  stress  intensity  factor 
equation  for  a finite  v/idth  center  cracked  plate,  and  then  apply  these  loads  to  the 
specimen.  Periodically  during  test,  crack  length  measurements  were  made  and  input  to 
the  computer  which  automatically  reduced  loads  in  accordance  with  the  stress  intensity 
factor  equation  defined  as 


where  P is  the  far-field  applied  load;  w and  B are  the  specimen  width  and  thickness, 
respectively;  and  a is  one-half  the  total  crack  length.  In  this  manner  testing  was  accom- 
plished under  quasi-constant  K conditions,  and  crack  length  data  were  updated  at  inter- 
vals necessary  to  prevent  desired  K values  varying  more  than  one  percent. 

In  addition  to  controlling  the  tests,  the  computers  were  also  used  to  store,  reduce 
and  output  data.  In  performing  a particular  test,  constant  amplitude  cycling  was  initial- 
ly performed  using  the  appropriate  R and  K2  values.  This  was  continued  while  collecting 
data  until  it  was  assured  that  a constant  growth  rate  had  been  reached.  The  desired  over- 
load and/or  underload  cycle  was  then  applied  followed  by  the  previously  applied  constant 
amplitude  condition.  During  this  latter  constant  amplitude  cycling,  crack  length  and 
cycle  data  were  collected,  examined  and  stored  until  it  was  assured  that  the  previous 
constant  amplitude  growth  rate  had  been  reached.  Except  for  cases  producing  a large 
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number  of  delay  cycles  and  hold  time  evaluations,  at  least  three  repeated  tests  were 
usually  performed.  For  these  repeated  tests,  crack  length  measurements  were  made  at 
the  same  cycle  increments  which  allowed  crack  length  data  to  be  averaged.  Upon  test 
completion,  the  stored  data  were  reduced  and  tabulated  by  the  computer  as  detailed  in 
Volume  II. 

For  approximately  50  percent  of  the  tests  producing  delay  cycles  in  excess  of 

5 . . 

10  , both  crack  tips  did  not  always  recover  at  the  same  rate  and  an  eccentric  crack  re- 
sulted. When  this  occurred,  testing  was  terminated  before  the  difference  in  finite  width 
correction  factors  for  the  symmetric  crack  and  eccentric  crack  invalidated  quasi-constant 
K conditions.  The  resulting  data  reflect  growth  history  of  the  one  crack  tip  that  had  at- 
tained constant  growth  rate.  All  other  data  represent  the  average  of  both  crack  tips. 

4.  RETARDATION  PARAMETER  DETERMINATION 

The  experimental  data  were  used  to  determine  the  number  of  delay  cycles,  Np, 
produced  by  the  overload  and/or  underload  cycle  and  the  associated  overload  affected 
zone  size,  a*,  defined  in  Figure  5.  Data  for  a selected  case  are  in  Figure  6 which  shows 
change  in  crack  length  after  overload  versus  number  of  cycles  after  overload  for  each  of 
three  repetitive  tests.  Also  shown  in  the  figure  are  average  data  for  the  three  tests. 
Figure  7 shows  incremental  Aa/ AN  versus  number  of  cycles  after  overload  derived  from 
the  basic  data  shown  in  Figure  6.  The  average  Aa/ AN  data  were  used  as  an  index  for 
determining  the  retardation  parameters,  and  for  the  case  illustrated  N^  was  56,000  cy- 
cles. Generally,  N^  could  be  determined  from  the  average  data  to  plus  or  minus  one 
data  interval  which  was  usually  either  1000  or  2000  cycles.  Figure  8 shows  Aa/  AN 
versus  change  in  crack  length  after  overload  and  an  a*  of  0.014  inch  corresponding  to 
the  56,000  delay  cycles  from  Figure  7.  Retardation  parameter  data  obtained  in  this 
manner  are  listed  in  Tables  V through  XVI  for  the  various  load  classes  and  stress  in- 
tensity factor  parameters  evaluated  during  the  experimental  program.  These  data  are 
illustrated  in  Section  IV  to  show  effects  of  the  various  load  classes  and  stress  intensity 
factor  parameters. 

5.  SHUT-OFF  RATIO  DETERMINATION 

In  performing  the  experimental  evaluations,  delay  cycle  and  overload  affected 
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zone  size  data  were  collected  at  0.5  Increasing  Increments  of  overload  ratio.  This 
was  continued  until  crack  arrest  occurred.  Additional  tests  were  then  performed  at 
overload  ratios  within  the  last  0.5  increment  to  zero-in  on  the  overload  shut-off  ratio, 
S , to  within  + O.I.  In  determining  overload  shut-off  ratio,  cycling  was  performed 

-9 

until  it  was  assured  that  crack  growth  rate  was  less  than  10  inch  per  cycle  which  cor- 

(2) 

responds  to  threshold  growth  rate  from  basic  da/dN  data''  . 
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SECTION  IV 

SUMMARY  OF  EXPERIMENTAL  DATA 
L EFFECTS  OF  OVERLOADS 

In  order  to  provide  baseline  data  for  determining  underload  effects,  several 
tests  were  performed  involving  tensile  overloads  without  underloads.  The  effect  of 
overload  on  crack  growth  history  is  illustrated  by  the  family  of  curves  in  Figure  9. 

These  curves  show  crack  extension  after  overload  versus  cycles  after  overload  for  dif- 
ferent values  of  overload  ratio  with  R constant  at  0.1  and  K2  constant  at  10  KSiyTrvT 
Delay  cycles  obtained  from  these  data  and  similar  data  for  R=0.3  and  0.5  are  summarized 
in  Figure  10.  Delay  cycles  increased  with  an  increase  in  overload  ratio,  S,  and  also 
Increased  as  increased  which  was  equivalent  to  an  increase  in  the  constant  amplitude 

stress  intensity  ratio,  R,  and  a decrease  in  AK.  The  overload  shut-off  ratio,  S , de- 
' so 

creased  as  AK  became  smaller  or  as  R became  larger.  In  this  figure,  and  similar  ones  to 
follow,  the  curves  have  been  drawn  to  project  toward  zero  delay  cycles  at  S=l  .0;  however, 
it  is  recognized  that  delay  may  not  exist  at  an  overload  ratio  greater  than  one.  The  over- 
load shut-off  ratio  points  were  plotted  on  the  10^  cycle  line  since  cycles  in  the  range  of 
5 X 10^  and  10^  were  applied  in  determining  shut-off  ratio.  Similar  data  showing  the  ef- 
fects of  K2  are  summarized  in  Figure  II  for  AK  constant  at  7 KSI\/inT  Delay  cycles  in- 
creased with  an  increase  in  overload  ratio,  but  decreased  as  K2  increased.  The  shut-off 
ratio  increased  as  K2  decreased. 

2.  EFFECTS  OF  UNDERLOADS 

A limited  number  of  tests  were  performed  to  determine  the  effect  of  a single 
underload  without  overload.  All  tests  were  performed  using  tension-tension  constant 
amplitude  cycles  with  tension,  zero  and  compression  underloads.  No  effect  on  crack 
growth  rate  was  observed  for  the  cases  evaluated.  Data  obtained  for  a high  compression 
underload  are  illustrated  in  Figure  12  which  shows  change  in  crack  length  after  underload 
versus  cycles  after  underload.  Figure  13  shows  corresponding  A a/  AN  versus  cycles  after 
underload.  Growth  rate  remained  essentially  constant  which  was  typical  for  all  under- 
load without  overload  tests. 

3.  EFFECTS  OF  UNDERLOAD  FOLLOWING  OVERLOAD 

The  effect  of  underload  on  decay  of  retardation  due  to  overload  is  illustrated 
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by  the  family  of  curves  in  Figure  14  which  shows  change  in  crack  length  and  cycles 
after  overload  for  K|,  and  held  constant  with  varying  underload.  Delay  cycles 
derived  from  such  data  are  shown  in  Figure  15  for  R constant  at  0.1  and  K2  constant  at 
10  KSI\ArT.  The  lower  curve  represents  no  underload.  A small  underload  produced  no 
marked  effect  on  delay  cycles  or  shut-off  ratio.  Underloading  into  compression  produced 
a relative  decrease  in  delay  cycles  and  an  increase  in  shut-off  ratio.  Note  also  that  no 
significant  underload  effect  was  shown  below  an  overload  ratio  of  two.  Similar  data 
for  R constant  at  0.3  are  shown  in  Figure  16.  There  was  a trend  that  tensile  and  zero 
underloads  decreased  delay  cycles  and  increased  shut-off  ratio^  but  it  could  only  be 
positively  concluded  for  the  compressive  underload.  As  was  observed  for  the  R=0.l  data, 
there  was  no  positive  underload  effect  at  S==2  and  below.  Data  for  R constant  at  0.5  are 
shown  in  Figure  17.  At  the  higher  R or  lower  AK,  the  effect  of  underloads  in  reducing 
delay  cycles  and  increasing  overload  shut-off  ratio  was  more  positive.  Also,  the  effect 
was  evident  at  lower  values  of  overload  ratio.  These  same  trends  were  observed  for 
K2='7.78  and  14  KSI\/in.  data  shown  in  F igures  18  and  19.  These  two  respective  figures 
represent  two  different  R values,  0.1  and  0.5,  but  a constant  AK  of  7 KSI\/Trv.  Data  for 
K^=0  from  these  figures  along  with  similar  data  for  K2=I0  are  shown  in  Figure  20  in  terms 
of  K|  and  delay  cycles.  Here  the  influence  of  an  increasing  K2  coupled  with  an  increasing 
^3  ~ ^4  shown  to  dramatically  decrease  delay  cycles  and  increase  shut-off  ratio. 

These  same  data  are  shown  in  Figure  21  in  terms  of  overload  ratio  to  demonstrate  the  abi- 
lity of  overload  ratio  to  normalize  the  K2  effect  at  constant  AK. 

In  Figures  18  and  19,  the  compressive  underloads  were  shown  to  be  approximate 
values  because  the  tension-compression  load  class  was  evaluated  in  terms  of  underload 
ratio,  consequently,  K^  varied  with  overload  ratio,  S.  It  was  found,  however, 

that  number  of  delay  cycles  was  relatively  insensitive  to  magnitude  of  K^.  This  is  il- 
lustrated in  Figure  22  which  shows  delay  cycle  variance  with  overload  and  underload 
ratios  for  K2=7.78  KSIv/irT.  S imilar  data  are  presented  in  Figure  23  for  K2=^I4  KSI\/Tn. 
Shut-off  ratio  data  were  not  obtained  for  K2=I4  KSIv^^Tru  since  Kj  reached  K critical 
before  shut-off  occurred.  The  same  general  results  were  obtained  for  tension-compression 
evaluations  in  terms  of  underload  ratio  at  K2=I0  KSI\/i^.  Also,  these  tension-compression 
evaluations  resulted  in  K|  values  greater  than  30  KSiyirT.  which  produced  stable  tear  and 
in  some  cases  severe  tunneling  during  the  overload  cycle.  This  is  subsequently  discussed 
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under  FRACTOGRAPHIC  EVALUATIONS.  Because  of  those  problems  evaluations  were 
performed  with  constant  at  -7.5  KSiyirT.  for  the  K2==I0  KSIv/in^.  tension-compression 
load  class. 

4.  EFFECTS  OF  UNDERLOAD  PRECEEDING  OVERLOAD 

In  addition  to  evaluating  underloads  following  overloads,  underloads  proceeding 
overloads  as  represented  by  the  zero-tension  and  compression-tension  load  classes  in 
Figure  I were  also  evaluated.  Figure  24  shows  the  compression-tension  data  along  with 
the  overload-underload  data  for  K2=I0  KSI\/Tr7.  and  R=0.l.  Delay  cycles  and  shut-off 
ratio  were  essentially  the  same  as  were  obtained  for  the  tension-tension  without  underload 
and  tension-zero  load  classes.  Only  the  tension-compression  load  class  showed  a signifi- 
cant effect.  Similar  data  for  R=0.3,  but  including  the  zero-tension  load  class,  are  shown 
in  Figure  25.  Again  no  definite  difference  was  observed  between  tension-tension,  tension- 
zero,  zero-tension  and  compression-tension.  A wider  variance  in  shut-off  ratio,  however, 
was  observed.  Figure  26  shows  data  for  R=0.5.  At  the  higher  R or  lower  AK;  the  under- 
load effect  of  the  tension-zero  load  class  became  evident,  but  a positive  difference  be- 
tween zero-tension,  compression-tension  and  tension-tension  without  underload  was  not 
present. 


Data  for  K2=7.78  KSiyin.  and  R=0.l  are  shown  in  Figure  27.  Again  the  tension- 
compression  effect  was  apparent,  and  tension-zero  showed  a reduction  in  delay  cycles 
above  S=2.5  and  an  increase  in  shut-off  ratio.  The  difference  between  tension-tension 
without  underload  and  compression  tension  was  not  positive.  Similar  data  for  K2=I4  KSI 
\/TnT  and  R=0.5  are  shown  in  Figure  28.  The  same  trends  were  observed,  but  underload 
effects  associated  with  tension-zero  were  evidenced  at  lower  overload  ratios.  From  all 
data,  it  was  evident  that  single  underloads  proceeding  overloads  produced  the  same  re- 
sults as  the  equivalent  tension-tension  case  without  underload. 

5.  OVERLOAD  SHUT-OFF  RATIO 

A summary  of  the  overload  shut-off  ratio  data  presented  in  the  previous  discus- 
sions is  presented  in  Figures  29  and  30.  Shut-off  ratio  variance  with  or  are  pre- 
sented in  Figure  29  for  tension-tension,  tension-zero,  and  tension-compression  load 
profiles.  S^^  is  shown  to  increase  as  the  underload  magnitude  increases.  It  is  also 

indicated  that  S increases  as  R decreases.  Shut-off  was  not  achieved  for  K«=I0 
so  2 
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and  14  KSI\/in.  when  followed  by  high  compressive  underloads  even  for  values  of  K| 
exceeding  the  stable  tear  threshold. 


Similar  data  are  shown  in  Figure  30  for  tension-tension  without  underload, 

zero-tension  and  compression-tension  load  profiles.  Underload  proceeding  overload 

produced  essentially  the  same  results  as  tension-tension  without  underload;  however, 

increased  as  R decreased.  In  both  Figures  29  and  30  there  was  no  distinguishable  effect 

(3) 

of  K2  variance  on  overload  shut-off  ratio.  This  supports  the  position  of  Hillberry^  ' that 
shut-off  ratio  is  independent  of  K2* 

6.  OVERLOAD  AFFECTED  ZONE  SIZE 

While  effects  of  the  different  load  classes  and  stress  intensity  parameters  on  delay 
cycles  were  of  primary  Interest,  the  overload  affected  zone  size,  a*,  was  also  of  interest 
since  the  models  require  an  a*  definition  for  delay  cycle  predictions.  It  was  previously 
shown  that  delay  cycles  Increased  with  an  increase  in  overload  ratio  for  any  particular 
set  of  conditions.  Also,  delay  cycles  increased  as  R increased  corresponding  to  a decrease 
in  AK.  The  overload  affected  zone  size  also  increased  as  overload  ratio  increased,  but 
decreased  as  R increased.  Figure  31  shows  overload  affected  zone  size  data  and  delay 
cycle  data  for  tension-tension  without  underload  and  K2  constant  at  10  KSIx/TirT  Both 
Np  and  a*  increased  as  overload  ratio  increased.  With  overload  ratio  constant,  delay 
cycles  increased  as  Increased,  but  a*  decreased.  Data  showing  these  same  trends  are 
shown  in  Figure  32  for  underloading  to  zero  following  overload,  and  in  Figure  33  for 
compressive  underload.  This  effect  on  overload  affected  zone  size  was  the  trend 
throughout  the  data  for  all  load  profile  classes  and  stress  intensity  factor  parameters. 

A positive  Influence  of  underload  magnitude  on  overload  affected  zone  size  could  not 
be  identified,  and  there  was  no  significant  difference  in  zone  size  data  for  underload 
proceeding  overload  and  tension-tension  without  underload.  Correlation  of  overload 
affected  zone  size  data  with  values  calculated  based  on  different  assumptions  are  con- 
tained in  a subsequent  section. 

7.  EFFECTS  OF  HOLD  PERIODS 

Evaluations  were  performed  to  determine  the  effects  of  hold  times  in  tension 
and  compression  on  delay  cycles.  Two  different  values  for  each  stress  intensity  fac- 
tor: K|,  K2f  Kg,  or  K^,  were  evaluated  for  hold  times  up  to  24  hours.  Values 
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for  each  stress  intensity  factor  were  selected  to  provide  delay  rather  than  shut-off 
based  on  data  for  no  hold  time.  The  family  of  growth  history  curves  in  Figure  34 
illustrates  the  increase  in  retardation  associated  with  increasing  hold  time  in  tension 
at  K| . Delay  cycle  data  obtained  for  hold  times  in  tension  are  summarized  in  Figures 
35  through  38.  In  every  case  a significant  increase  in  delay  cycles  was  observed  as 
hold  time  at  K|  was  increased.  There  was  no  indication  that  stable  conditions  had  been 
reached  at  the  24  hours  hold  period;  consequently,  a further  increase  in  delay  cycles 
could  be  expected  for  longer  hold  periods.  The  data  in  Figure  38  for  K^-5  KSl\/in • 
indicate  that  the  overload  shut-off  ratio  was  being  approached  at  the  24  hour  hold  time. 
Delay  cycles  for  K|=25  KSIv/TrT.  in  Figure  35  and  for  K2=7.78  KSiyTrT.  in  Figure  36  were 
also  approaching  the  shut-off  range  at  24  hours.  It  would  be  reasonable  to  assume  that 
shut-off  would  have  been  reached  for  these  cases  with  a short  extension  in  hold  time. 
Approximately  paralleling  results  were  obtained  for  each  value  of  the  same  stress  intensity 
factor.  The  paralleling  results  in  Figure  38  for  K^=5  KSI\/in~.  or  overload  without  under- 
load and  for  K^=l  KSI\/irK  representing  underloading  in  tension  indicated  that  underloading 
after  holding  in  tension  had  no  effect.  Although  delay  cycles  increased  significantly  with 
hold  time,  there  was  no  distinguishable  effect  on  overload  affected  zone  size  as  evidenced 
by  data  in  Table  XV. 


Data  for  hold  time  in  compression  are  illustrated  by  the  family  of  crack  length 
versus  cycle  curves  shown  in  Figure  39.  The  data  indicate  a small  decrease  in  delay 
cycles,  and  that  stable  conditions  were  essentially  reached  within  the  first  hour. 

These  results  were  typical  for  all  cases  evaluated  as  shown  by  the  delay  cycle  versus 
hold  time  data  in  Figures  40  through  43. 


8.  FRACTOGRAPHIC  EVALUATIONS 

When  evaluating  the  tension-compression  load  class  for  constant  values  of  U / 

c 

K|/K^,  it  was  necessary  to  exceed  a K|  of  30  KSiyTn".  Even  then  the  overload  shut-off 
ratio  was  not  obtained  for  many  cases  because  K|  was  incremented  upward  until  it  reached 
K critical  during  the  overload  cycle.  As  K|  was  increased  above  30  KSiyTr7.  relatively 
large  plastic  zones  were  first  observed  at  the  crack  tip  and  branching  or  deviation  from 
the  original  crack  path  during  recovery  was  also  observed.  These  conditions  became  more 
severe  when  K|  exceeded  35  KSI\/in^,  and  in  many  cases  overload  affected  zone  sizes 
determined  from  the  crack  growth  data  were  considerably  larger  than  the  specimen  thickness. 
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The  fracture  surfaces  of  several  specimens  from  the  tension-compression  load  class 
were  microscopically  examined  to  classify  the  crack  front  geometry  for  several  values 
of  K|  • The  photographs  in  Figure  44  show  fracture  surfaces  containing  overloads  of  25 
and  30  KSIv/TnT.  At  lOX  magnificat  ion  the  25  KSIx/irT.  overload  produced  arrest  marks 
on  the  fracture  surface  which  were  hardly  detectable.  The  30  KSI\/Tr7.  mark  was  more 
pronounced  and  both  crack  fronts  exhibited  a slight  curvature. 

In  Figure  45  similar  photographs  are  shown  for  K|=35  and  40  KSI\/Tn^  At  the  35 
KSIv/irT.  overload^  the  mark  was  even  more  pronounced.  The  40  KSI\/Tri.  overload,  how- 
ever, caused  considerable  crack  extension  and  the  onset  of  tunneling  was  evident.  Sizeable 
shear  lips  were  also  observed  at  the  specimen  surfaces.  The  tunneling  and  shear  lips  are  very 
evident  for  the  K|=49  KSI\/in^  case  shown  in  Figure  46.  The  specimen  surface  is  also  shown 
in  this  figure  to  illustrate  the  deviation  of  crack  path  observed.  The  lines  on  the  specimen 
surface  represent  the  previously  described  grill  against  which  crack  growth  measurements 
were  made.  Based  on  crack  growth  rate  measurements  the  case  illustrated  in  Figure  46 
produced  an  overload  affected  zone  size  of  0.490  inch  or  approximately  twice  the  specimen 
thickness. 

Based  on  the  fracture  surface  observations,  delay  cycle  and  overload  affected  zone 
size  data  reported  for  Kj  values  greater  than  35  KSiyirv.  are  considered  invalid  and  of  ques- 
tionable value  from  a linear  elastic  fracture  mechanics  point  of  view.  The  questionable 
cases  are  limited  to  data  in  Tables  XI 1 and  XIII  where  K|  exceeded  35  KSI  \/in • Rather  than 

rationalize  data  for  these  cases,  additional  data  were  collected  with  constant  at  -7.5 
KSiyTnT.  in  order  to  evaluate  the  tension-compression  load  class.  These  data  are  reported 
in  Table  XIV  and  were  contained  in  several  of  the  previous  data  illustrations. 
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SECTION  V 

ANALYTICAL  METHODS 


1 . SUMMARY 

Analytical  predictions  of  delay  cycles  were  made  for  the  same  cases  evalua- 
ted experimentally.  In  making  these  predictions  four  basic  elements  were  required. 

They  were  (I)  baseline  crack  growth  rate  data  or  crack  growth  rate  equation,  (2) 
stress  intensity  factor  solution,  (3)  applied  load  spectrum,  and  (4)  the  load-interaction 
models  and  their  assumed  overload  affected  zones.  The  accuracy  of  the  fatigue  crack 
growth  prediction  depends  upon  the  accuracy  of  all  individual  elements.  A good  crack 
growth  retardation  model  alone  does  not  necessarily  guarantee  a good  prediction. 

Since  we  have  high  confidence  in  the  stress  intensity  factor  solution.  Equation  (I), 
for  a center  cracked  specimen;  and  high  precision  of  load  spectrum  control,  these 
two  parameters  offer  minimal  uncertainty  in  the  predictions.  The  other  two  parameters, 
namely  the  baseline  crack  growth  equation  and  the  crack  growth  retardation  models  and 
their  assumed  overload  affected  zones,  play  very  important  roles  in  prediction  of  delay 
cycles.  In  this  section,  the  crack  growth  rate  equations,  retardation  model  formulations 
and  overload  affected  zones  used  in  subsequent  delay  cycle  predictions  are  discussed. 

2.  CRACK  GROWTH  RATE  EQUATIONS 

Three  different  crack  growth  rate  equations  from  the  literature  were  employed 
in  making  delay  cycle  predictions  with  each  of  the  three  retardation  models.  The 
equations  were  independently  developed  for  22I9-T85I  aluminum  alloy  and  were  the 
Hall^^^,  Grumman  and  Walker  equations.  Each  represented  slightly  different 
mathematical  formulations  fitted  to  basic  da/dN  versus  AK  data.  The  specific  mathe- 
matical formulations  employed  were  as  follows: 

(I)  Hall  Crack  Growth  Rate  Equation 

Crack  growth  rate  data  generated  using  three  stress  ratios  (R=0.1, 

0.3  and  0.5)  were  fitted  using  the  following  equation: 

W'  (2) 
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Values  of  c , and  n^,  were  defermined  by  least  square 
ri  n ri 

fitting  the  data  by  assuming  K^|^“l  *5  KSI  The  computed 

values  were 

c = 0.34  X 10  ^ 

H 

m = 0 . 84 

"h 

(2)  Grumman  Crack  Growth  Rate  Equation 

The  data  generated  for  various  R.  ratios  were  fitted  to  a modi- 
fied Elber  equation  of  the  form: 

- c[(l  + qR)4K]" 


dN 


(3) 


A least  squares  procedure  was  used  to  fit  the  data  to  the  above 

equation,  and  the  best  curve  fit  produced  the  following  constants: 

c = 1 .96  X 10  ^ 

n = 3.34 

q = 0.6 

R = 0.5 
c 

where  R is  the  stress  ratio  cutoff  value,  i.e.,  if  R>R  , R=R  . 
c c c 

(3)  Walker  Crack  Growth  Rate  Equation 

The  Walker  equation  has  the  following  general  form: 

\A/  L 


dN  W max ' 

When  the  data  were  fitted  to  this  equation  the  following 
constants  were  obtained: 


(4) 


c - I *72  X 10 

W 

= 0.3014 
W 

nw=  3.415 


-9 


The  unit  for  the  stress  intensity  factor  is  KSI>/irv!,  and  the 
crack  growth  rate  is  inch/cycle. 


In  order  to  understand  the  trends  those  equations  predicted,  the  da/dN  vs. 
AK  computed  using  the  three  equations  were  plotted  in  Figures  47  to  49  for  R=0.l, 
0.3  and  0.5,  respectively.  As  seen  from  those  figures,  for  AK  > 3 KSI\/in .,  the 
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difference  in  growth  rate  between  the  Hall  and  Grumman  equations  was  negligible 
for  all  three  R values.  But,  when  4K  is  less  than  3 KSI  \/ i n . , such  difference  be 
comes  more  significant.  The  smaller  the  4K,  the  larger  the  difference  in  growth 
rate,  and  the  difference  decreases  as  R increases.  The  primary  reason  for  the  dif- 
ference is  that  the  Hall  equation  is  anchored  by  a threshold  value.  The  Walker 
equation  matched  the  Grumman  equation  very  well  for  R— 0,1;  however,  when  R 
increases,  the  growth  rate  predicted  by  the  Walker  equation  becomes  much  higher 
than  that  of  the  Grumman.  However,  when  R ratio  is  less  than  0.1,  the  rate  compu- 
ted using  the  Hall  equation  becomes  larger  than  those  computed  using  the  Grumman 
and  Welker  equations.  Figure  50  shows  the  relationship  between  da/dN  and  4K 
for  R=  -0.3.  As  seen  from  this  figure,  when  4K  is  greater  than  3 KSI  \J in .,  the 
growth  rate  computed  using  the  Hall  equation  is  higher  than  those  computed  using 
the  Grumman  and  Walker  equations. 

3.  RETARDATION  MODELS 

It  is  well  known  that  normal  crack  growth  rate  under  constant  amplitude 
loading  changes  if  the  load  application  is  preceded  by  a loading  of  different  ampli- 
tude. A tensile  overload  produces  permanent  plastic  deformation  at  the  crack  tip 
which  delays  the  crack  growth  of  subsequent  low  load  cycles.  The  beneficial  effect 
of  tensile  overloads  in  extending  crack  life,  however,  may  be  altered  if  overloads 
are  coupled  with  underloads.  Current  analytical  models  deal  with  the  decay  of  this 
beneficial  effect  as  a function  of  crack  tip  progress  through  an  "overload  affected 
zone"  . In  order  to  evaluate  the  ability  of  current  models  to  predict  underload  effects, 
delay  cycle  predictions  were  made  for  the  same  cases  evaluated  experimentally.  Pre- 
dictions were  made  using  three  different  retardation  models:  (I)  generalized  Wheeler 
model (2)  generalized  Willenborg  model and  (3)  Grumman  closure  model 
Formulations  of  the  three  aforementioned  crack  growth  retardation  models  are  pre- 
sented and  discussed  in  Appendix  A. 

4.  OVERLOAD  AFFECTED  ZONE 

The  overload  affected  zone,  a*,  due  to  the  application  of  the  overload  is 
defined  as  the  distance  ahead  of  the  crack  tip,  over  which  the  reduction  in  fatigue 
crack  growth  rate  from  that  of  constant  amplitude  loading  is  observed.  To  date. 
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fhis  zone  has  been  characterized  by  the  plastic  zone  size  created  by  the  overload  cycle. 
The  Wheeler  model  relates  a*  to  the  plane  strain  plastic  zone  radius  while  the  Willen- 
borg  and  closure  models  estimate  a*  as  a function  of  the  plane  stress  plastic  zone  radius. 

/o\  /p\ 

Some  other  investigators,  e.g.,  Probst^  , Von  Euw^  , utilized  variations  of  the  overload 
created  plastic  zone  diameter  (depending  upon  the  fraction  of  the  cross  section  exhibiting 
shear  lips)  to  achieve  good  correlation  with  measured  overload  affected  zone. 


All  three  retardation  models  predict  a return  to  steady  state  (constant  amplitude) 
crack  growth  rate  subsequent  to  the  application  of  the  overload  as  soon  as  the  total 
crack  growth  increment  since  the  application  of  the  overload  equals  the  overload  af- 
fected zone,  i.e.,  4a=a*.  If  the  condition  is  expressed  in  terms  of  the  stress  intensity 

factor,  tfien  the  generalized  Wheeler  and  Willenborg  models  predict  the  termination  of 

* * 

retardation  effect  due  to  overload  when  K and  K are  equal . K is  the  ef- 

max  max  max 

fective  stress  intensity  factor  after  overload  and  K is  the  stress  intensity  factor  of  the 
' max 

constant  amplitude  maximum.  The  generalized  closure  model  will  predict  the  same  thing 
when  K^2  equal . being  closure  stress  intensity  factor  for  constant  ampli- 

tude and  K is  the  effective  closure  stress  intensity  factor  after  overload.  The  overload 

^ "k 

affected  crack  length  can  be  determined  from  the  condition  K =K  which  results 

max  max 

in  the  following  equation' 


K 2 

f max  ^ 

max 


For  the  generalized  closure  model 


where  Zp.. , defined  earlier  as  the  overload  created  load  interaction  zone, 

. (9) 

has  the  following  general  form' 


2 

Zql 
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(8) 


The  factor  a can  be  expressed  in  the  following  general  form 

I 


(3). 


a - 


2X  + 4/T  (1-X) 


Y 


7T 


where  X represents  the  fraction  of  the  cross  section  exhibiting  shear  lips, 
and  Y equals  I or  2 depending  upon  use  of  the  overload  plastic  zone  radius 
or  diameter. 


Figures  51,  52  and  53  show  a comparison  between  the  measured  overload  af- 
fected crack  length,  and  the  computed  a*  using  Equation  (5)  and  various  factors  for 
the  different  overload  classes.  The  factor  ce  used  in  computing  a*  ranged  from  0.0563 
representing  the  plane  strain  plastic  zone  radius  to  0.3183  representing  the  plane  stress 
plastic  zone  diameter.  As  seen  from  these  figures,  the  majority  of  data  indicates  that 
the  measured  overload  affected  zone  decreases  with  the  increase  of  the  stress  intensity 
factor  ratio,  R.  This  suggests  that  the  a*  may  depend  not  only  on  the  maxima  of  the 
overload  and  constant  amplitude  stress  intensity  factors  (K|  and  K2)/  but  also  the  minimum 

of  the  constant  amplitude  stress  intensity  factor,  K^. 

o 
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SECTION  VI 
MODEL  EVALUATIONS 


L COMPUTER  PROGRAM 

A computer  program  was  prepared  for  the  purpose  of  predicting  delay  cycles 
due  to  the  application  of  a single  overload  and/or  underload.  The  delay  cycles  to 
be  predicted  depend  not  only  upon  the  model,  but  also  upon  the  crack  growth  rate 
equation  and  the  overload  affected  zone.  For  fairness  of  model  evaluations,  the 
program  was  prepared  to  predict  delay  cycles  for  each  of  the  three  models  combined 
with  each  of  the  three  crack  growth  equations.  Also,  the  program  allowed  predictions 
using  the  plastic  zone  radius  and  diameter  for  both  plane  strain  and  plane  stress  condi- 
tions. 


A numerical  integration  approach  similar  to  the  one  used  by  Gallagher  and 
Hughes^^  was  employed  to  estimate  the  number  of  applied  load  cycles  required  to 
propagate  the  crack  from  the  position  right  after  the  application  of  overload  to  a posi- 
tion one  overload  affected  zone  ahead  of  the  crack  tip.  Four  different  increments, 
ranging  from  1/25  to  1/250  of  a*  were  used  in  the  numerical  integration  to  study  the 
rate  of  convergence  of  integration,  and  results  are  tabulated  in  Table  XVII  for  a case 
where  K|=I5,  units  in  KSI\Ain.  The  delay  cycles  were 

computed  using  the  Grumman  growth  rate  equation  and  the  plane  stress  plastic  zone 
radius  as  the  overload  affected  zone.  As  can  be  seen  from  the  table,  the  rate  of  con- 
vergence was  very  rapid.  Since  the  computer  time  required  to  execute  the  program 
was  small,  100  integration  steps  were  used  in  the  subsequent  predictions. 

2.  MODEL  SENSITIVITY  TO  SHUT-OFF  RATIO 

The  only  parameter  capable  to  account  for  the  underload  effect  in  the  mathe- 
matic formulation  of  the  generalized  Wheeler  and  Willenborg  models  is  the  overload 
shut-off  ratio,  S^^.  For  a given  constant  amplitude  stress  intensity  factor  ratio,  ex- 
perimental results  indicate  that  the  overload  shut-off  ratio  is  dependent  upon  the 
magnitude  of  underload  or  compressive  load.  The  larger  the  magnitude  of  compressive 
load,  the  higher  the  overload  shut-off  ratio. 

For  a high  compressive  load  immediately  following  the  tensile  overload,  there 
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may  not  be  a shut-off  in  crack  growth.  Even  if  there  is  a shut-off,  the  shut-off 
overload  may  be  much  higher  than  stable  tear  threshold  value  or  the  net  section 
stress  (in  the  plane  containing  the  crack  surface)  may  exceed  the  yield  strength  of  the 
material . In  this  case  the  value  of  shut-off  ratio  may  not  be  too  meaningful  from  linear 
elastic  fracture  mechanics  point  of  view. 

Since  the  overload  shut-off  ratio  is  required  in  the  generalized  Wheeler  and 
Willenborg  models  and  generation  of  such  ratios  for  all  levels  of  underloads  is  impractical, 
a study  was  conducted  to  determine  the  effect  of  the  overload  shut-off  ratio  on  predicted 
delay  cycles.  It  was  found  that  the  predicted  delay  cycles  is  fairly  insensitive  to  the 
overload  shut-off  ratio  if  the  overload  cycle  minimum  is  close  to  the  constant  amplitude 
cycle  minimum.  However,  for  a tension-compression  overload,  such  effect  is  more  pro- 
nounced. Figure  54  shows  the  predicted  delay  cycles  normalized  by  the  measured  delay 
cycles  as  a function  of  assumed  overload  shut-off  ratios  for  K|=20,  K2=I0,  K^=  -K|  and 
various  K^.  As  can  be  seen  from  the  figure,  the  effect  of  on  the  predicted  delay 
cycles  is  more  significant  in  the  generalized  Wheeler  model  than  the  generalized  Willen- 
borg model . This  effect  remains  fairly  constant  for  various  constant  amplitude  stress 
intensity  factor  ratios,  R.  The  important  thing  to  be  pointed  out  is  that  if  the  lower  over- 
load shut-off  ratio  is  used  in  the  prediction,  the  delay  cycles  may  be  over  estimated.  For 
the  purpose  of  predicting  delay  cycles  due  to  overload  and/or  underload,  the  curves  repre- 
senting the  overload  shut-off  ratios  as  shown  in  Figures  29  and  30  were  used.  These  were 
obtained  by  fitting  the  curves  through  data  generated  for  each  R ratio. 

The  generalized  closure  model  does  not  depend  upon  the  overload  shut-off  ratio. 

It  accounts  for  the  underload  or  compressive  load  effect  through  the  closure  factor,  C^, 
of  the  overload  cycle. 

3.  MODEL  SENSITIVITY  TO  THRESHOLD  STRESS  INTENSITY  RANGE 

Other  than  the  overload  shut-off  ratio,  both  the  generalized  Wheeler  and 
generalized  Willenborg  models  also  require  the  value  of  threshold  stress  intensity  fac- 
tor range,  Since  experimental  scatter  has  been  observed  in  data  for  determining 

a study  was  conducted  to  determine  the  sensitivity  of  the  models  to  this  parameter. 
The  predicted  delay  cycles  using  the  generalized  Wheeler  and  generalized  Willenborg 
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models  are  presented  as  a function  of  threshold  AK  in  Figure  55  for  both  tension- 
tension  and  tension-compression  overloads.  The  compressive  load  chosen  was  the' 
one-half  of  the  tension  overload.  Two  levels  of  overload  (S=l.5  and  3,0)  and  the 
Grumman  crack  growth  rate  equation  were  used  in  the  prediction.  In  each  overload 
case,  the  prediction  was  made  by  holding  the  overload  shut-off  ratio  constant  at  mea- 
sured value  and  varying  threshold  AK  from  I .!  and  1.9.  As  seen  from  this  figure,  the 
effect  of  threshold  AK  on  predicted  delay  cycles  increases  with  the  increase  of  overload 
ratio.  The  generalized  Willenborg  model  exhibited  a low  sensitivity  to  the  threshold 
value  while  the  generalized  Wheeler  model  is  more  sensitive  to  the  threshold  value. 

The  average  threshold  K obtained  in  Reference  2 is  about  1.5  for  R=0  and  was  used  in 

max 

the  analytical  predictions  of  delay  cycles.  For  applied  load  ratio,  R,  different  from  zero; 
threshold  AK  was  set  equal  to  1.5  (l-R). 

4.  EVALUATION  OF  THE  GENERALIZED  WHEELER  MODEL 

As  discussed  in  Section  V,  the  predicted  delay  cycles  due  to  the  application  of 
overload  does  not  depend  only  on  the  model  and  the  crack  growth  rate  equation,  but  also 
on  the  overload  affected  zone,  a*.  For  any  given  constant  amplitude  and  overload  stress 
intensity  factors,  the  a*  is  directly  proportional  to  the  load  interaction  zone  created  by 
the  overload,  • Both  the  original  and  generalized  Wheeler  model  assumed  that  the 
Zql  was  equal  to  the  overload  created  plane  strain  plastic  zone  radius.  Although  most  of 
the  measured  overload  affected  crack  lengths  were  larger  than  the  corresponding  ones  com- 
puted using  Equation  (5)  by  assigning  z^j^  as  the  overload  created  plane  strain  plastic  zone 
radius,  these  computed  a**s  were  used  in  the  prediction  of  delay  cycles  due  to  various 
types  of  overloads. 

The  correlations  of  measured  and  predicted  delay  cycles  using  the  generalized 
Wheeler  model  and  the  three  different  crack  growth  rate  equations  are  presented  in 
Figures  56  through  61  for  22I9-T85I  aluminum  plates  subjected  to  various  types  of  over- 
load and/or  underload.  In  all  figures,  an  arrow  connected  to  the  symbol  Indicates  that 
the  predicted  delay  cycles  is  outside  the  range  of  the  graphic  scale.  The  solid  line  in 
all  these  figures  represents  perfect  correlation,  and  the  broken  lines  represent  a scatter 
factor  of  two.  The  data  used  In  correlations  Include  three  different  constant  amplitude 
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stress  intensity  factor  ratios  (R=0.l,  0.3  and  0.5)  and  three  different  constant  ampli- 
tude stress  intensity  maxima  (K2=10,  7.78  and  14  KSI  \/Tru).  The  overload  ratio,  S, 
ranges  from  1 .5  to  3.0  for  all  load  classes  except  tension-compression  in  which  S ranged 
from  1.5  to  4.5.  In  general  the  higher  the  overload  ratio,  the  larger  the  delay  cycles. 

Figure  56  shows  the  correlations  of  measured  and  predicted  delay  cycles  due 
to  tension-tension  overloads.  When  the  generalized  Wheeler  model  is  used,  the  pre- 
dicted delay  cycles  using  the  Hall  growth  rate  equation  is  smaller  than  the  corresponding 
one  predicted  using  the  Grumman  or  Walker  equations.  Except  for  a smaller  overload 
ratio  (S=l  .5),  the  predicted  delay  cycles  is  larger  than  the  actual  measured  delay  cycles, 
especially  for  a high  overload  ratio  even  though  a smaller  than  measured  overload  af- 
fected zone  was  used  in  predictions.  If  a larger  overload  affected  zone  (as  the  one  mea- 
sured) is  used,  then  the  predicted  delay  cycles  will  be  even  higher  than  the  current 
predictions.  Figure  57  shows  similar  correlations  of  delay  cycles  due  to  tension-zero 
overloads  for  K2=7.78,  10  and  14  KSIv/lrT.  Figures  58  and  59  present  the  delay  cycles 
due  to  the  application  of  tension-compression  overload  cycles  for  K2=I0  and  K2"=7.78 
and  14  KSI  \An .,  respectively.  By  comparing  the  results  shown  in  Figures  57  through  59 
with  that  of  56,  one  can  see  that  the  presence  of  underload  or  compressive  load  imme- 
diately after  the  overload,  does  not  alter  the  trend  of  the  prediction  on  delay  cycles. 
When  the  sequence  of  the  overload  cycle  is  reversed,  i.e.  the  overload  cycle  is  applied 
in  an  order  of  zero-tension  or  compressive-tension,  the  retarded  crack  growth  behavior 
observed  is  very  similar  to  those  of  tension-tension  overloads.  The  generalized  Wieeler 
model  does  not  account  for  the  sequence  of  overload  cycle.  Nevertheless,  the  generali- 
zed Wheeler  model  was  used  to  predict  the  delay  cycles  due  to  zero-tension  and  compres- 
sion-tension overloads,  and  the  results  are  presented  in  Figures  60  and  61 . As  was  shown 
for  the  other  overload  classes,  the  generalized  Wheeler  model  overestimates  the  delay 
cycles  for  all  cases  where  the  overload  ratio  is  higher  than  1 .5. 

In  general,  the  predicted  delay  cycles  by  the  generalized  Wheeler  model  together 
with  the  usage  of  the  overload  created  plane  strain  plastic  zone  radius  in  computing  the 
overload  affected  zone,  is  significantly  higher  than  the  measured  delay  cycles  whenever 
the  overload  stress  intensity  factor  exceeds  1 .5  times  the  constant  amplitude  stress  in- 
tensity factor  maximum.  At  an  overload  ratio  of  1 .5,  most  of  the  predicted  delay  cycles 
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are  within  a factor  of  two  of  the  actual  measured  values. 


5.  EVALUATION  OF  THE  GENERALIZED  WILLENBORG  MODEL 

As  mentioned  earlier,  the  generalized  Willenborg  model  assumes  that  the 
overload  created  interaction  zone,  equal  to  the  overload  created  plane  stress 

plastic  zone  radius.  The  computed  overload  affected  zone,  a*,  using  the  overload 
created  plane  stress  plastic  zone  radius  correlates  better  with  the  data  than  the  one 
computed  using  the  overload  created  plane  strain  plastic  zone  radius  as  noted  in  Figures 
51  through  53.  For  the  purpose  of  evaluating  the  generalized  Willenborg  model,  the 
overload  affected  crack  length,  a*,  computed  using  the  overload  created  plane  stress 
plastic  zone  radius  together  with  three  crack  growth  rate  equations  discussed  in  Section 
V were  used  in  predicting  the  delay  cycles  due  to  various  types  of  overloads.  Based 
upon  the  current  study,  the  predicted  delay  cycles,  using  the  generalized  Willenborg 
model  and  the  computed  a*  using  the  overload  created  plane  stress  plastic  zone  radius 
is  significantly  lower  than  the  measured  delay  cycles.  In  most  cases,  the  predictions  are 
too  conservative.  A better  correlation  was  observed  when  the  plane  stress  plastic  zone 
diameter  was  used  to  compute  the  overload  affected  crack  length.  Therefore,  the  corre- 
lations of  predicted  and  measured  delay  cycles  using  the  generalized  Willenborg  model 
are  presented  for  both  the  a*  computed  using  the  overload  created  plane  stress  plastic  zone 
radius  and  diameter.  Figures  62  through  67  show  correlations  of  measured  delay  cycles 
and  delay  cycles  predicted  using  as  the  plane  stress  plastic  zone  radius.  In  these 
figures,  the  solid  line  represents  perfect  correlation,  and  the  broken  lines  represent  a 
scatter  factor  of  two.  Examination  of  these  figures  show  that  predictions  were  generally 
unconservative  for  all  load  profiles. 

Similar  correlations  are  shown  in  Figures  68  through  73  using  os  the  plane 
stress  plastic  zone  diameter.  Figure  68  shows  the  correlations  of  predicted  and  measured 
delay  cycles  using  the  generalized  Willenborg  model  and  three  different  crack  growth 
rate  equations  for  22I9-T85I  aluminum  alloy  plate  subjected  to  tension-tension  overloads. 
Similar  results  on  delay  cycles  due  to  tension-zero  overloads  are  shown  in  Figure  69. 

As  seen  from  these  two  figures  compared  with  Figures  62  and  63,  a better  correlation 
will  be  achieved  if  the  overload  affected  zone  is  approximated  by  the  overload  created 
plane  stress  plastic  zone  diameter.  All  predictions  are  within  a factor  of  two  of  the 
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measured  values,  regardless  of  which  crack  growth  rate  equation  is  employed. 

Figure  70  shows  correlations  of  measured  and  computed  delay  cycles  due  to 
tension-compression  overloads.  Although  the  correlation  is  not  as  good  as  that  of 
tension-tension  and  tension-zero  overload  cases,  the  overall  correlation  is  still  good. 

As  discussed  earlier  in  this  section,  one  of  the  main  problems  in  predicting  the  delay 
cycles  due  to  tension-compression  overloads  is  the  selection  of  a proper  overload 
shut-off  ratios.  The  overload  shut-off  ratios  shown  in  Figures  29  and  30  were  the  best 
estimations  based  upon  the  available  data.  For  high  compressive  load,  the  shut-off 
overload  is  much  higher  than  stable  tear  threshold  value.  In  order  to  enhance  the 
evaluation  of  the  model  in  predicting  delay  cycles  due  to  tens  ion- compressive  over- 
loads, additional  data  including  the  overload  shut-off  ratio  were  generated  for  K2=I0 
and  Kg=-7.5  KSl/u^  The  correlations  of  measured  and  predicted  delay  cycles  using 
the  generalized  Willenborg  model  for  these  additional  cases  is  shown  in  Figure  71. 

Again,  good  correlation  was  obtained  if  the  a*  was  computed  using  the  plane  stress 
plastic  zone  diameter.  When  the  application  of  underload  or  compressive  load  precedes 
the  overload,  the  observed  overload  shut-off  ratio  does  not  change  significantly  with  the 
increase  in  magnitude  of  compressive  load  as  seen  in  Figure  30.  Therefore,  one  could 
anticipate  that  the  retarded  crack  growth  behavior  due  to  compression-tension  overload 
might  be  close  to  that  of  tension-tension  overload.  This  is  confirmed  by  the  results  pre- 
sented in  Figures  72  and  73  which  show  the  correlations  of  measured  and  predicted  delay 
cycles  due  to  zero-tension  and  compressive-tension  overloads,  respectively. 

In  general,  if  the  overload  created  load  Interaction  zone,  is  assumed  equal 

to  the  overload  created  plane  stress  plastic  zone  radius  as  proposed  by  the  Willenborg  model, 
then  the  delay  cycles  predicted  by  the  generalized  Willenborg  is  almost  always  less  than 
the  actual  measured  delay  cycles,  except  for  some  of  the  tension-compression  overload 
cases.  The  crack  growth  rate  equations,  proposed  by  Hall,  Bell  and  Walker  for  22I9-T85I 
aluminum  alloy  plate,  do  not  show  a consistent  trend  in  the  predicted  delay  cycles  due 
to  various  overloads  and/or  underloads.  It  is  almost  impossible  to  rank  one  equation  to 
the  others.  In  general,  any  one  of  these  three  equations  may  be  used  alone  to  evaluate 
the  models. 

When  the  overload  affected  zone  was  computed  using  the  overload  created  plane 
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stress  plastic  zone  diameter,  almost  all  the  predicted  delay  cycles  using  the  generalized 
Wlllenborg  model  were  within  a factor  of  two  of  the  actual  measured  delay  cycles. 

6.  EVALUATION  OF  THE  CLOSURE  MODEL 

The  closure  model  assumes  that  the  overload  affected  zone,  a*,  Is  equal  to  the 
overload  created  plane  stress  plastic  zone  radius.  Unlike  the  generalized  Wheeler  and 
generalized  Wlllenborg  models,  the  closure  model  does  not  depend  upon  the  overload 
shut-off  ratio  to  account  for  the  underload  or  compressive  load  effect.  Such  effect  was 
accounted  for  through  the  closure  factor,  Cp  from  which  the  closure  stress  intensity 
factor,  K^,  was  determined.  Upon  application  of  the  overload,  if  the  computed 
is  larger  than  the  constant  amplitude  stress  intensity  factor  maximum,  then  the  closure 
model  predicts  that  the  crack  arrest  occurs. 

The  correlations  of  measured  and  predicted  delay  cycles  using  the  closure  model 
and  various  crack  growth  rate  equations  for  22I9-T85I  aluminum  alloy  plates  subjected  to 
various  types  of  overloads  and/or  underloads  are  shown  in  Figures  74  through  79.  Figure 
74  presents  delay  cycle  correlations  due  to  tension-tension  overloads.  As  seen  from  this 
figure,  the  predicted  delay  cycles  using  the  Hal!  equation  are  always  less  than  the  corres- 
ponding ones  predicted  by  the  Grumman  and  Walker  equations.  Predictions  using  the 
Grumman  and  Walker  equations  show  almost  no  difference.  For  an  overload  ratio,  S,  less 
than  or  equal  to  2.0,  the  predicted  delay  cycles  using  the  closure  model,  and  the  Hall 
crack  growth  rate  equation,  are  within  a factor  of  two  of  the  measured  delay  cycles. 
However,  for  S>2,0,  the  predicted  delay  cycles  are  much  higher  than  the  measured 
values,  and  crack  arrest  is  predicted  regardless  of  which  crock  growth  rate  equation  is  used 
whenever  Si>2,5.  Similar  results  were  observed  for  tension-zero  overload  cases  which  are 
shown  in  Figure  75.  Figures  76  and  77  show  the  correlations  of  measured  and  predicted 
delay  cycles  due  to  tension-compression  overloads  for  ^2”^^  KSiyinT  and  K2~7.78  and 
14  KSIn/w^.,  respectively.  Figures  78  and  79  present  similar  results  on  delay  cycles  due 
to  zero-tension  and  compression-tension  overloads,  respectively.  Except  for  compression- 
tension  overload  cases,  the  trend  of  delay  cycle  predictions  by  the  closure  model  is  the 
same  as  that  for  tension-tension  overload  cases.  The  predicted  delay  cycles  due  to 
compression-tension  overloads  are  less  than  that  of  measured  values. 

In  general,  when  the  closure  model  (together  with  the  assumption  that  the  over- 
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load  effected  zone  equals  the  overload  created  plane  stress  plastic  zone  radius)  is 
used  to  predict  the  delay  cycles  due  to  any  type  of  overload  and/or  underload,  the 
best  results  will  be  achieved  if  the  Hall  crack  growth  rate  equation  is  employed.  If 
the  overload  ratio  is  less  than  or  equal  to  two,  the  predicted  delay  cycles  are  within  a 
factor  of  two  of  the  measured  delay  cycles.  If  the  Grumman  or  Walker  crack  growth 
rate  equation  is  used,  together  with  the  closure  model,  to  predict  the  delay  cycles  due 
to  overloads,  the  predicted  delay  cycles  will  consistently  be  higher  than  the  measured 
values.  When  the  overload  stress  intensity  factor  exceeds  2,5  times  the  constant  ampli- 
tude stress  intensity  factor  maximum,  the  closure  model  will  predict  that  the  crack  arrest 
occurs,  while  there  are  finite  delay  cycles  (less  than  half  million  cycles)  measured. 

7.  COMPARISON  OF  MODELS 

As  discussed  during  the  evaluation  of  each  model,  the  predicted  delay  cycles 
due  to  overload  and/or  underload  is  dependent  upon  the  overload  affected  zone  size 
chosen  in  the  prediction.  It  was  found  that  the  best  results  will  be  achieved  by  pairing 
the  model  and  the  overload  affected  zone  as  follows: 

1)  the  generalized  Wheeler  model  and  the  plane  strain  plastic  zone 
radius, 

2)  the  generalized  Willenborg  model  and  the  plane  stress  plastic 
zone  diameter, 

3)  the  closure  model  and  the  plane  stress  plastic  zone  radius 

In  order  to  have  a direct  comparison  between  models,  the  predictions  using  the 
Hall  crack  growth  rate  equation  and  the  above  pairings  of  the  model  and  overload  affect- 
ed zone,  are  tabulated  in  Tables  XVIII  through  XXVII  for  various  types  of  overload  and/or 
underload  conditions.  Those  results  are  also  presented  in  Figures  80  through  86.  The 
reason  the  Hall  equation  was  chosen  is  that  it  correlates  with  available  data  better  than 
the  Grumman  and  Walker  equations  in  the  region  where  AK  is  less  than  3 KSi  v/irT. 

When  AK  is  higher  than  3 KSI  >/in/,  the  difference  in  growth  rate  predicted  by  three 
equations  is  practically  negligible. 

Figure  80  shows  the  correlations  of  measured  and  predicted  delay  cycles  using 
the  Hall  crack  growth  rate  equation  and  three  aforementioned  retardation  models  for 
22I9-T85I  aluminum  alloy  plates  subjected  to  tension-tension  overloads.  The  data  in- 
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eludes  three  constant  amplitude  stress  intensity  factor  maxima  {K^7 10  and 
14  KSIv/^irT.)  and  three  applied  load  ratios  (R=0J,  0.3  and  0.5).  The  overload  ratio 
ranges  from  1.5  to  3.0.  Within  the  same  overload  class,  the  higher  the  overload 
ratio,  the  larger  the  delay  cycles.  As  seen  from  Figure  80  when  the  overload  ratio 
is  greater  than  2.5,  the  closure  model  predicts  the  occurrence  of  crack  arrest  (or 
infinite  delay  cycles),  while  the  generalized  Wheeler  model  predicts  delay  cycles 
greater  than  two  times  the  measured  value  in  most  cases,  even  though  a smaller  than 
measured  a*  was  used  in  the  prediction.  The  generalized  Willenborg  model  gives  the 
best  prediction  for  S^^2.5.  When  the  overload  ratio  is  below  2.0,  both  the  generalized 
Willenborg  model  and  the  closure  mode!  give  a prediction  within  a factor  of  two  of 
measured  delay  cycles.  In  this  range,  the  prediction  using  the  generalized  Wheeler 
model  is  not  as  good  as  the  other  two  models.  Similar  correlations  of  delay  cycles  due 
to  tension~zero  overloads  are  presented  in  Figure  81.  Those  results  show  practically  no 
difference  from  those  of  tension-tension  overload  class. 

Correlations  of  measured  and  predicted  delay  cycles  due  to  tension-compression 
overloads  are  shown  in  Figure  82  for  K2“I0  and  in  Figure  75  for  K2=7.78  and  14.  As 
discussed  earlier  in  this  section,  both  the  generalized  Wheeler  and  generalized  Willenborg 
models  require  the  overload  shut-off  ratio  for  each  overload  condition  while  the  closure 
model  does  not  depend  upon  such  ratio.  Experimental  determination  of  accurate  overload 
shut-off  ratios  for  tension-compression  overload  class  is  impractical . The  overload  shut- 
off ratios  used  to  obtain  the  predicted  delay  cycles  shown  in  Figures  82  and  83  were  ex- 
trapolated values.  As  seen  from  Figures  82  and  83,  for  tension-compression,  the  correla- 
tion using  the  generalized  Wheeler  and  generalized  Willenborg  models  are  not  as  good  as 
the  ones  obtained  for  tension-tension  and  tension-zero  overload  classes.  For  an  overload 
ratio  less  than  or  equal  to  2.0,  the  predictions  using  the  closure  model  correlate  fairly 
good  with  the  measured  values.  However,  when  the  overload  ratio  is  greater  than  or 
equal  to  2.5,  the  prediction  using  the  closure  model  is  completely  no  good.  Addi- 
tional data  were  generated  for  the  compressive  underload  constant  at  -7,5  KSIvAnT, 
and  results  of  measured  and  predicted  delay  cycles  using  the  three  retardation  models  are 
shown  in  Figure  84.  With  availability  of  the  overload  shut-off  ratio  for  -7.5  KSI\/irK, 
the  predictions  using  the  generalized  Willenborg  model  excellently  correlate  with  the 
measured  values,  while  the  closure  model  gives  a compatible  good  prediction  with  the 
generalized  Willenborg  model  for  S ^ 2 ,0,  but  bad  predictions  whenever  S >2  .5. 
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Correlations  of  measured  and  predicted  delay  cycles  due  to  zero-tension  and 
compression-tension  overloads  are  shown  in  Figures  85  and  86,  respectively.  As  was 
obtained  for  the  other  overload  classes,  the  best  prediction  was  obtained  by  using  the 
generalized  Willenborg  model . Both  the  generalized  Wheeler  model  and  the  closure 
model  give  a fairly  good  prediction  for  lower  overload  ratios  but  bad  predictions  for 
high  overload  ratios. 

To  summarize  the  model  comparison,  one  may  conclude  that: 

1)  When  overload  ratio  is  greater  than  or  equal  to  2.5  the  delay 
cycles  predicted  by  the  closure  model  is  completely  no  good 
regardless  of  what  type  of  overload;  the  delay  cycles  predicted 
by  the  generalized  Wheeler  model  is  much  larger  than  the 
measured  value,  while  the  generalized  Willenborg  model  gives 
the  best  prediction  which  is  within  a factor  of  two  of  measured 
delay  cycles. 

2)  When  overload  ratio  is  less  than  or  equal  to  2.0,  all  three 
models  provide  a reasonable  prediction  except  for  the 
tension-compression  overload  class  having  high  compressive 
load.  For  this  overload  case,  the  closure  model  gives  the 
best  prediction. 
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SECTION  VII 

GENERAL  CONCLUSIONS 


1)  Underloads  following  overloads  decreased  the  retardation  effect  produced  by 
the  overload  and  increased  the  overload  shut-off  ratio.  These  effects  became  more 
apparent  as  the  constant  amplitude  minimum  stress  intensity  factor  was  increased  while 
holding  other  parameters  constant.  For  high  compressive  underloads  following  over- 
loads, a physically  meaningful  overload  shut-off  ratio  may  not  exist. 

2)  Underloads  proceeding  overloads  produced  approximately  the  same  number  of 
delay  cycles,  overload  affected  zone  size  and  overload  shut-off  ratio  as  the  equivalent 
tension-tension  case  with  the  underload  truncated. 

3)  In  this  investigation  the  constant  amplitude  maximum  stress  intensity  factor 
was  varied  over  an  approximate  range  of  7 to  14  KSI  \/ In . Data  obtained  within  that 
range  indicate  that  the  overload  shut-off  ratio  is  independent  of  the  constant  amplitude 
maximum . 

4)  Hold  periods  in  tension  increased  delay  cycles,  and  stabilization  was  not 
apparent  after  24  hours  which  was  the  longest  period  evaluated.  Although  delay  cy- 
cles Increased  with  hold  time  there  was  no  measurable  effect  on  overload  affected 
zone  size.  It  was  also  indicated  that  holding  in  tension  decreases  the  overload  shut- 
off ratio. 

5)  Hold  periods  in  compression  resulted  in  a small  decrease  in  delay  cycles, 
and  the  effect  was  essentially  stable  after  one  hour. 

6)  Based  upon  the  experimental  data  generated  under  this  program,  it  was  found 
that  the  overload  affected  zone,  a*,  depends  not  only  on  the  maxima  of  overload  and 
constant  amplitude  load  cycles,  but  also  on  the  minimum  of  the  constant  amplitude  cycle. 

7)  For  overload  ratios  higher  than  two  times  the  constant  amplitude  stress  intensity 
factor  maximum,  the  Grumman  closure  model,  in  general,  gives  a very  unconservative 
prediction.  In  most  cases,  the  closure  model  predicts  that  crack  arrest  occurs  while 
finite  delay  cycles  were  measured. 
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8)  The  generalized  Wheeler  model  consistently  predicts  higher  delay  cycles  than 
the  actual  measured  delay  cycles,  even  though  a much  smaller  than  measured  overload 
affected  zone  was  used  in  the  prediction.  This  implies  that  the  effective  AK  computed 
by  the  model  is  much  smaller  than  the  actual  value. 

9)  Excellent  correlation  of  measured  and  predicted  delay  cycles  can  be  achieved 
by  using  the  generalized  Willenborg  model,  provided  that  a larger  than  the  measured 
overload  affected  zone  is  used  in  the  prediction.  Physically,  it  implies  that  the  effec- 
tive AK  computed  by  the  generalized  Willenborg  model  is  larger  than  the  actual  value. 
This  is  just  opposite  to  that  computed  by  the  generalized  Wheeler  model . 

10)  The  predicted  delay  cycles  due  to  a single  overload  is  highly  dependent  upon 
the  size  of  overload  affected  zone  chosen.  If  the  subsequent  overload  is  applied  before 
the  crack  tip  extends  outside  the  retarded  region  created  by  the  first  overload,  then  the 
retarded  growth  rate  predicted  should  be  less  sensitive  to  the  a*.  In  general,  such 
sensitivity  will  decrease  with  the  decrease  of  period  between  overloads, 

11)  The  effect  of  three  different  crack  growth  rate  equations  discussed  earlier  on 
the  predicted  delay  cycles  was  small.  This  was  partially  due  to  the  high  constant  ampli- 
tude stress  intensity  factor  ranges  ( 4K=5,  7,  9 KSI\/in^)  used  in  the  program.  When 

AK  was  greater  than  3 KSI\/in.,  the  growth  rate  computed  using  the  three  different 
equations  show  practically  no  difference.  If  a lower  AK  is  used  in  the  analysis,  sub- 
stantial difference  in  the  predicted  delay  cycles  by  using  different  growth  rate  equation 
may  be  anticipated. 

12)  The  overload  shut-off  ratio  is  fairly  sensitive  to  the  magnitude  of  compressive 
load  following  the  tensile  overload.  If  the  generalized  Wheeler  model  or  generalized 
Willenborg  model  is  used  in  the  life  prediction  of  fatigue  crack  growth  under  varying 
amplitude  spectrum  load,  a choice  of  proper  overload  shut-off  ratios  may  be  a problem, 
especially  for  a new  material  where  the  shut-off  ratio  data  are  not  available. 
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SECTION  VIII 
RECOMMENDATIONS 


Based  on  results  of  this  program,  areas  have  been  identified  as  recommended 
future  efforts  directed  toward  improving  crack  growth  prediction  methodology.  These 
are  briefly  discussed  below  for  consideration  in  future  programs  sponsored  by  the  Air 
Force  or  industry. 

I)  The  overload  affected  zone  size  formulation  plays  a very  important  role  in 

crack  growth  predictions.  Generally,  current  formulations  consider  only  the  maxima 
of  the  overload  and  constant  amplitude  cycles  as  controlling  parameters.  The  ma[ority 
of  data  reported  herein  indicate  that,  with  all  other  parameters  constant,  the  measured 
overload  affected  zone  size  decreases  with  an  increase  in  stress  intensity  factor  ratio,  R. 
This  suggests  that  a*  may  depend  not  only  on  the  maxima  of  the  overload  and  constant 
amplitude  stress  intensity  factors,  but  also  the  minimum  of  the  constant  amplitude  cycle. 

In  order  to  improve  current  prediction  methodology  or  to  develop  a better  load-interaction 
model,  development  of  a mathematical  model  to  characterize  the  overload  affected  zone 
more  accurately  is  necessary.  Some  additional  data  are  considered  necessary  to  identify 
the  influence  of  key  parameters  and  guide  the  formulation  of  such  a model , 

2)  Current  retardation  models  do  not  account  for  hold  time  effects  on  subsequent 
crack  growth,  but  current  data  indicate  that  they  should  be  modified  to  do  so.  Accord- 
ing to  current  data,  hold  time  in  tension  increases  the  delay  cycles  but  does  not  practically 
change  the  overload  affected  zone  size.  Incorporation  of  the  overload  shut-off  ratio  as 

a function  of  hold  time  in  the  generalized  Wheeler  and  generalized  Willenborg  models 
or  the  closure  stress  intensity  factor  as  a function  of  hold  time  in  the  closure  model  appears 
to  be  the  proper  direction  for  hold  time  accountability  in  tension.  Holding  in  compres- 
sion resulted  in  a small  decrease  in  delay  cycles,  and  additional  data  are  required  to 
determine  the  need  or  an  approach  for  accountability.  Also,  an  extension  of  the  data 
reported  herein  on  holding  in  tension  will  be  required  to  guide  model  modifications. 

3)  Since  both  the  generalized  Wheeler  and  generalized  Willenborg  models  require 

overload  shut-off  ratios,  it  is  recommended  that  procedures  for  determining  S and 

so 

mathematically  representing  in  the  models  be  established.  Data  reported  herein 


32 


indicate  that  underload  preceeding  overload  can  be  adequately  represented  by  the 
equivalent  overload  case  v/ith  the  underload  truncated,  and  that  is  independent 
of  the  constant  amplitude  maximum.  Consequently,  the  important  parameters  for  con- 
sideration are  limited  to  the  magnitude  of  underload  following  overload  and  the  constant 
amplitude  load  ratio.  Mathematically  representing  may  be  more  straight  forward  for 
the  generalized  Willenborg  model  since  it  is  less  sensitive  to  variance  than  the 
generalized  Wheeler  model . Consequently,  it  may  be  possible  to  adequate!  y represent 

S as  a function  only  of  underload  in  the  generalized  Willenborg  model . 
so  ' 

4)  All  data  generated  under  this  program  were  for  single  overload/underload  cases. 
Data  for  periodic  loading  or  simple  blocked  spectrum  should  be  developed  as  the  next 
step  in  model  evaluations. 

5)  In  crack  growth  analysis,  a good  load- interaction  model  alone  does  not  neces- 
sarily guarantee  a good  prediction.  The  prediction  is  also  highly  dependent  upon  the 
baseline  crack  growth  rate  data  and  stress  intensity  factor  solution.  Efforts  in  developing 
a better  model  or  improving  the  existing  models  should  couple  with  the  development  of 
more  reliable  baseline  data  as  well  as  method  for  presenting  these  data. 
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FIG.  1(d)  TENSION-COMPRESSION  FIG.  1(e)  COMPRESSION-TENSION 


Figure  I.  Load  Profile  Classifications 
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S=K|/K2/  Overload  Ratio 

R ~ Constant  Amplitude  Stress  Intensity  Ratio 

U“K|/K^,  Underload  Stress  Intensity  Ratio 
U^=K|/K^,  Compressive  Stress  Intensity  Ratio 

AK=K2"K2,  Constant  Amplitude  Stress  Intensity  Range 
Figure  2.  Stress  Intensity  Factor  Parameter  Definitions 
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Figure  5.  Stress  Intensity  Factor  Parameters 
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CYCLES  AFTER  OVERLOAD 

Figure  9.  Effect  of  Overload  Ratio  on  Crack  Growth  for  Tension-Tension 


Load  Closs  Without  Underload,  K«=I0,  R=0.l,  K„=K. 

2 3 4 


<3=3 


Figure  10.  Effect  of  Overload  Ratio,  S,  and  Constant  AmpI itude  Stress 

Intensity  Ratio,  R^K^/K-,  on  Delay  Cycles  for  Tension-Tension 

Load  Class  Without  Underload,  K«=I0,  K^*K. 

2 3 4 


DELAY  CYCLES 

Figure  II.  Effect  of  Overload  Ratio,  S,  ond  on  Delay  Cycles  for 

Constant  Amplitude  Stress  Intensity  Range,  AK=K-  -K^,  of  7 
Without  Underloads 


THOUSANDS  OF  CYCLES  AFTER  UNDERLOAD 


Figure  12,  Crack  Growth  History  For  Underload  Without  Overload 
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THOUSANDS  OF  CYCLES  AFTER  UNDERLOAD 
igure  13.  Crack  Growth  Rate  History  For  Underload  Without  Overload  - From  Figure  12  Data 


DELAY  CYCLES 


Figure  16.  Effect  of  Overload  Ratio  and  Underload  on  Delay  Cycles  for 
R=0.3  and  K2=10 


DELAY  CYCLES 

Figure  17,  Effect  of  Overload  Ratio  and  Underload  on  Delay  Cycles  for 
R=0.5  and  K„=I0 
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Figure  19.  Effect  of  Overload  Ratio  and  Underload  on  Delay  Cycles  for  R=0.5  and  <^=14 
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Figure  24.  Effect  of  Load  Class  on  Deloy  Cycles  for  Constont  Amplitude 
Stress  Intensity  Rotio,  R=K^K^.  of  0,1  and  K^^IO 


Figure  25.  Effect  of  Lood  Class  on  Deloy  Cycles  for  Constant  Amplitude 
Stress  Intensity  Ratio,  R=K2/K2  of  0.3  and  K2=I0 


Figure  26,  Effect  of  Load  Class  on  Delay  Cycles  for  Constant  Amplitude 
Stress  Intensity  Ratio,  R=K^/K2,  of  0.5  and  K2=I0 
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CYCLES  AFTER  OVERLOAD,  THOUSANDS 
Figure  34.  Effect  of  Hold  Time  in  Tension  on  Crock  Growth  History 


HOLD  TIME  (i  Kj,  HOURS 


Figure  35.  Effect  of  K|  and  Hold  Time  at  K|  on  Delay  Cycles 
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HOLD  TIME  AT  K|,  HOURS 

Figure  36.  Effect  of  and  Hold  Time  at  K|  on  Delay  Cycles 
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HOLD  TIME  AT  HOURS 

Figure  37.  Effect  of  and  Hold  Time  at  K|  on  Delay  Cycles 
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gure  38.  Effect  of  and  Hold  Time  at  K|  on  Delay  Cycles 
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Figure  41 , Effects  of  and  Hold  Time  at  on  Delay  Cycles 
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Figure  43.  Effects  of  and  Hold  Time  at  on  Delay  Cycles 
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Figure  44.  Fracture  Surface  Features  Produced  by  Tensile  Overloads  of 
25  and  30  KSiyin.,  Magnification  = 10 
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Tensile  Overload 
Location  - K.=35  KSI  yiru 


Figure  45.  Fracture  Surface  Features  Produced  by  Tensile  Overloads  of 
35  and  40  KSI  Magnification  - 10 
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Figure  46.  Fracture  Surface  and  Specimen  Surface  Showing  Features  Produced  by  a 
Tensile  Overload  of  49  KSI\/ln. 
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da/dN,  INCH/CYCLE 


Figure  50.  da/dN  vs.  4K  Calculated  for  R=  -0.3  Using  Three  Different  Growth 
Rate  Equations  for  22I9-T85I  Aluminum  Alloy 
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MEASURED  OVERLOAD  AFFECTED  ZONE,  INCH 


MEASURED  OVERLOAD  AFFECTED  ZONE,  INCH 


PLANE  STRESS  DIAMETER 
PLANE  STRESS  RADIUS 

plane  strain  diameter 

PLANE  STRAIN  RADIUS 


PREDICTED  OVERLOAD  AFFECTED  ZONE,  INCH 


Figure  52  . Correlations  of  Predicted  and  Measured  Overload 

Affected  Zones  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Tension-Compression  Overloads 


MEASURED  OVERLOAD  AFFECTED  ZONE,  INCH, 


PREDICTED  DELAY  CYCLES/MEASURED  DELAY  CYCLES 


OVERLOAD  SHUT-OFF  RATIO,  S 

so  _ 

Figure  54.  Effect  of  Overload  Shut-Off  Ratio  on  Predicted  Delay  Cycles 
for  K|=20,  K2=I0,  K^=20  and  Various 
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PREDICTED  DELAY  CYCLES 


MEASURED  DELAY  CYCLES 
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PREDICTED  DELAY  CYCLES 


Figure  56.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Wheeler  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  SubjeOted  to 
Tension-Tension  Overloads 
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PREDICTED  DELAY  CYCLES 


Figure  57.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Wheeler  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Tension-Zero  Overloads 
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Figure  58.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Wheeler  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Tension-Compression  Overloads  and  K2=I0  KSI  sf\n. 


69 


O hall  equation 

A GRUMMAN  EQUATION 
□ WALKER  EQUATION 
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PREDICTED  DELAY  CYCLES 


Figure  60.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Wheeler  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subject  to 
Zero-Tension  Overloads 
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Figure  61.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Wheeler  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Compression-Tension  Overloads 
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MEASURED  DELAY  CYCLES 


Figure  63.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using 

the  Generalized  Willenborg  Model  and  Various  Crack  Growth 
Rate  Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected 
to  Tension-Zero  Overloads 
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MEASURED  DELAY  CYCLES 
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Figure  64.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using 

the  Generalized  Willenborg  Model  and  Various  Crack  Growth 
Rate  Equations  for2219-T851  Aluminum  Alloy  Plate  Subjected 
to  Tension-Compression  Overloads 
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Figure  65.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Willenborg  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Tension-Compression  Overloads  and  K2=I0  and  K^=  -7.5  KSI\/inT 
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MEASURED  DELAY  CYCLES 
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Figure  66.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Willenborg  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Zero-Tension  Overloads 
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Figure  69.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using 

the  Generalized  Willenborg  Model  and  Various  Crack  Growth 
Rate  Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected 
to  Tension-Zero  Overloads 
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Figure  70.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using 

the  Generalized  Willenborg  Model  and  Various  Crack  Growth 
Rate  Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected 
to  Tension-Compression  Overloads 
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Figure  71.  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Willenborg  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Tension-Compression  Overloads  and  K«=I0  and  Kj-=  -7.5  KSI\/in, 
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MEASURED  DELAY  CYCLES 


Figure  73,  Correlations  of  Predicted  and  Measured  Delay  Cycles  Using  the 
Generalized  Willenborg  Model  and  Various  Crack  Growth  Rate 
Equations  for  22I9-T85I  Aluminum  Alloy  Plate  Subjected  to 
Compression-Tension  Overloads 
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Figure  74  • Correlations  of  Measured  and  Predicted  Delay  Cycles 
Using  the  Closure  model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Tension-Tension  Overloads 
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MEASURED  DELAY  CYCLES 


PREDICTED  DELAY  CYCLES 


Figure  75  . Correlations  of  Measured  and  Predicted  Delay  Cycles 
Using  the  Closure  Model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Tension- Zero  Overloads 
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MEASURED  DELAY  CYCLES 
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PREDICTED  DELAY  CYCLES 


Figure  76  , Correlations  of  Measured  and  Predicted  Delay  Cycles 
Using  the  Closure  Model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Tension-Compression  Overloads  and 

K2»  10  Ksiyin 
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MEASURED  DELAY  CYCLES 


10^  10^  10^  10^ 

PREDICTED  DELAY  CYCLES 

Figure  77  . Correlations  of  Measured  and  Predicted  Delay  Cycles 


Using  the  Closure  Model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Tension-Compression  Overloads  and  • 14 
and  7.78  Ksi  v/lrT.  / 
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MEASURED  DELAY  CYCLES 
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Figure  78  . Correlations  of  Measured  and  Predicted  Delay  Cycles 
Using  the  Closure  Model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Zero-Tension  Overloads 
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MEASURED  DELAY  CYCLES 


PREDICTED  DELAY  CYCLES 

Figure  79  . Correlations  of  Measured  and  Predicted  Delay  Cycles 
Using  the  Closure  Model  and  Various  Crack  Growth 
Rate  Equations  for  2219-T851  Aluminum  Alloy  Plates 
Subjected  to  Compression-Tension  Overloads 
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MEASURED  DELAY  CYCLES 


Figure  80  • Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Tension-Tension  Overloads 
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MEASURED  DELAY  CYCLES 


PREDICTED  DELAY  CYCLES 


Figure  81  . Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Tension- Zero  Overloads 
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Figure  82  . CorrelaMons  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Tension-Compression  Overloads  and 
K„=  10  Ksi/liT 
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Figure  83  . Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Tension-Compression  Overloads  and 
= 7.78  and  14  Ksi,/kT 
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Figure  84  • Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Tension-Compression  Overloads  and 
K„=  10  and  K.  = -7.5  Ksi/l^ 
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Figure  85  . Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for  2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Zero-Tension  Overloads 
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Figure  86  . Correlations  of  Predicted  and  Measured  Delay  Cycles 

Using  the  Hall  Crack  Growth  Rate  Equation  and  Various 
Retardation  Models  for '2219-T851  Aluminum  Alloy 
Plates  Subjected  to  Compression-Tension  Overloads 
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TABLE  V 

RETARDATION  PARAMETER  DATA  FOR  22I9-T85I  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 
DESICCATED  AIR  FOR  TENSION-TENSION  LOAD  CLASS 
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(1)  Overload  shut-off  ratio 

(2)  Data  for  one  crack  tip. 
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RETARDATION  PARAMETER  DATA  FOR22I9-T85I  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 

DESICCATED  AIR  FOR  TENSION-ZERO  LOAD  CLASS 
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(1)  Data  for  one  crack  tip. 

(2)  Overload  shut-off  ratio. 
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(I)  Overload  shut-off  ratio. 


TABLE  XI 
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TABLE  XII 


RETARDATION  PARAMETER  DATA  FOR  2219-T85I  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 
DESICCATED  AIR  FOR  TENSION-COMPRESSION  LOAD  CLASS 
><2=10 


Specimen 

No. 

S 

R 

— 

AK 

U 

c 

•<1 

•<2 

s 

^D 

a*,  Inch 

2-L-7 

1.0 

0.1 

9 

-.67 

10 

10 

1 

-15 

0 

0 

2-L-7 

1.0 

0.5 

5 

-.67 

10 

10 

5 

-15 

0 

0 

4-L-I4 

1.5 

0.1 

9 

-1.0 

15 

10 

1 

-15 

4500 

0.017 

4-L-I4 

1.5 

0.1 

9 

-2.0 

15 

10 

1 

-7.5 

4500 

0.016 

6-L-3 

1.5 

0.3 

7 

-1.0 

15 

10 

3 

-15 

5500 

0.015 

6-L-3 

1.5 

0.3 

7 

-2.0 

15 

10 

3 

-7.5 

6500 

0.017 

2-L-I2 

1.5 

0.5 

5 

-1.0 

15 

10 

5 

-15 

6000 

0.007 

5-L-2 

1.5 

0.5 

5 

-2.0 

15 

10 

5 

-7.5 

6000 

0.007 

2-L-7 

2.0 

0.1 

9 

-1.0 

20 

10 

1 

-20 

14000 

0.048 

2-L-7 

2.0 

0.1 

9 

-2.0 

20 

10 

1 

-10 

16000 

0.050 

3-L-3 

2.0 

0.3 

7 

-1.0 

20 

10 

3 

-20 

12000 

0.029 

5-L-I7 

2.0 

0.3 

7 

-2.0 

20 

10 

3 

-10 

14000 

0.026 

l-L-19 

2.0 

0.5 

5 

-1.0 

20 

10 

5 

-20 

12000 

0.012 

3-L-I8 

2.0 

0.5 

5 

-2.0 

20 

10 

5 

-10 

16000 

0.013 

4-L-I4 

2.5 

0.1 

9 

-1.0 

25 

10 

1 

-25 

26000 

0.051 

6-L-I3 

2.5 

0.1 

9 

-2.0 

25 

10 

1 

-12.5 

29000 

0.056 

l-L-19 

2.5 

0.3 

7 

-1.0 

25 

10 

3 

-25 

28000 

0.057 

6-L-12 

2.5 

0.3 

7 

-2.0 

25 

10 

'3 

-12.5 

26000 

0.050 

l-L-19 

2.5 

0.5 

5 

-1.0 

25 

10 

5 I 

-25 

24000 

0.015 

6-L-12 

2.5 

0.5 

5 

-2.0 

25 

10 

5 

-12.5 

37000 

0.023 

4-L-ll 

3.0 

0.1 

9 

-1.0 

30 

10 

1 

-30 

51000 

0.103 

2-L-I2 

3.0 

0.1 

9 

-2.0 

30 

10 

1 

-15 

83000 

0.091 

4-L-ll 

3.0 

0.3 

7 

-1.0 

30 

10  1 

3 i 

-30 

66000 

0.077 

5-L-7 

3.0 

0.3 

7 

-2.0 

30 

10 

3 

-15 

88000 

0.090 

5-L-7 

3.0 

0.5 

5 

-1.0 

30 

10 

5 

-30 

60000 

0.041 

5-L-I4 

3.0 

0.5 

5 

-2.0 

30 

10 

5 

-15 

II 1000 

0.043 

3-L-7 

3.5 

0.1 

9 

-1.0 

35 

10 

1 

-35 

116000 

0.174 

3-L-7 

3.5 

0.1 

9 

-2.0 

35 

10 

1 

-17.5 

251000 

0.127  (1) 

6-L-I5 

3.5 

0.3 

7 

-1.0 

35 

10 

3 

-35 

11 1000 

0.169 

2-L-3 

3.5 

0.3 

7 

-2.0 

35 

10 

3 

-17.5 

272000 

0.211  (1) 

6-L-I5 

3.5 

0.5 

5 

-1.0 

35 

10 

5 

-35 

147000 

0.096  (1) 

2-L-3 

3.5 

0.5 

5 

-2.0 

35 

10 

5 

-17.5 

635000 

0.261 

2-L-I7 

4.0 

0.1 

9 

-1.0 

40 

10 

1 

-40 

95000 

0.206 

l-L-5 

4.0 

0.1 

9 

-2.0 

40 

10 

1 

-20 

187000 

0.321 

5-L-2I 

4.0 

0.3 

7 

-1.0 

40 

10 

3 

-40 

165000 

0.224  (!) 

5-L-2I 

4.0 

0.3 

7 

-2.0 

40 

10 

3 

-20 

707000 

0.319  (1) 

5-L-2! 

4.0 

0.5 

5 

-1.0 

40 

10 

5 

-40 

326000 

0.240 

2-L-6 

4.0 

0.5 

5 

-2.0 

40 

10 

5 

-20 

1430000 

0.303  (!) 

(I)  Data  for  one  crack  tip. 
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TABLE  XIII 

RETARDATION  PARAMETER  DATA  FOR  22I9-T85I  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 
DESICCATED  AIR  FOR  TENSION-COMPRESSION  LOAD  CLASS 

K=\4  and  K^=7.78 


0.032 

0.027 

0.018 

0.034 

0.063 

0.092 

0.110 

0.180  (1) 
0.533  (1) 

0.490 
0.703 

0.012 
0.012 
0.016 
0.015 
0.034 
0.024 
0.049 
0.047 
0.037 
0.108  (1) 
0.061  (1) 
0.109  (1) 

0.107  (1) 

D 

z 

8500 

7500 

6500 

14000 

28000 

45000 

50000 

151000 

246000 

276000 

316000 

8000 

9000 

14000 

16000 

36000 

40000 

42000 

92000 

102000 

172000 

266000 

246000 

452000 

(2) 

(2) 

LO 

-21 

-10.5 

-28 

-14 

-35 

-17.5 

-42 

-21 

-49 

-24.5 

-56 

-11.67 

-5.84 

-15.56 

-7.78 

-19.45 

-9.72 

-23.34 

-II  .67 

-27.23 

-13.62 

-31 .12 

-15.56 

-35.01 

-15.95 

-38.90 

CO 

00  00  00  00  00  00  00  00  00  00  00  00  00  00  00 
KhsKfxrvKKKKKKKKKrs 

KKhvKKr\KI\KI\r\  ooooooooodooooo 

CN 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

7.78 

7.78 

7.78 

7.78 

7.78 

7.78 

7.78 
7.78 
7.78 
7.78 
7.78 
7.78 
7.78 

i 7.78 
1 7.78 

21 

21 

28 

28 

35 

35 

42 

42 

49 

49 

56 

11.67 

11.67 

15.56 

15.56 

19.45 

19.45 

23.34 

23.34 

27.23 

27.23 

31.12 

31 .12 

35.01 

31.90 

38.90 

U 

c 

— CNI  — CSI  — CN  — CN—  CN  — — CN  — CN  — CN  — CN  CN  CN  CN  

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 t 1 t 1 

iNlNINrNlNlNlNlNhvrNrN  rNlNrNrNKlNlNlNrNlNlNrNrNlNlN 

oc 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

O.I 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

I 

' 0.1 

1 0.1 

1.5 

1.5 
2.0 
2.0 

2.5 

2.5 
3.0 

3.0 

3.5 

3.5 

4.0 

1 .5 

1.5 

2.0 
2.0 

2.5 

2.5 
3.0 

3.0 

3.5 

3.5 

4.0 

4.0 

4.5 

4.1 
5.0 

c 

<D 

,§  0 
u 2 
0) 

CL 

LO 

6-L-20 

6-L-20 

3-L-I4 

2- L-9 

3- L-I4 
3-L-I7 
6-L-l 

5- L-3 

6- L-l 

1- L-5 

2- L-I7 

6-L-2I 

6-L-2I 

5-L-20 

5- L-20 

3- L-I5 

3- L-I5 
2-L-I9 

4- L-2I 
2-L-I9 
2-L-I6 
2-L-I9 

6- L-l 
6-L-l 
4-L-7 
4-L-7 
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(1)  Data  for  one  crack  tip 

(2)  Overload  shut-off  rati 
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(1)  Data  from  basic  program. 

(2)  Overload  shut-off  ratio. 


TABLE  XV 


RETARDATION  PARAMETER  DATA  FOR  22I9-T851  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 
DESICCATED  AIR  FOR  TENSION-TENSION  LOAD  CLASS 


HOLD  TIME  EFFECTS 


Specimen 

No. 

s 

R 

U 

V 

Time  K., 

Hours 

% 

a*,  Inch 

6-L-I4 

2.0 

0.1 

9 

20 

20 

10 

1 

1 

0 

16000 

0.042 

6-L-I4 

2.0 

0.1 

9 

20 

20 

10 

! 

1 

0.0014 

16000 

0.033 

6-L-I4 

2.0 

0.1 

9 

20 

20 

10 

1 

1 

0.25 

22000 

0.042 

6-L-I4 

2.0 

0.1 

9 

20 

20 

10 

1 

1 

1.0 

22000 

0,035 

5-L-M 

2.0 

0.1 

9 

20 

20 

10 

1 

1 

24 

28000 

0.038 

6-L-I4 

2.5 

0.1 

0 

25 

25 

10 

1 

0 

43000 

0.043 

6-L-I4 

2.5 

0.1 

9 

25 

25 

10 

1 

1 

0.25 

70000 

0.044 

6-L-I4 

2.5 

0.1 

9 

25 

25 

10 

1 

1.0 

89000 

0.029 

6-L-I4 

2.5 

O.I 

9 

25 

25 

10 

1 

1 

4.0 

102000 

0.048 

(1) 

2-L-6 

2.5 

0.1 

9 

25 

25 

10 

1 

24 

127000 

0.032 

2-L-4 

2.57 

0.13 

6.78 

20 

25.7 

7.78 

1 

0 

72000 

0.027 

2-L-4 

2.57 

0.13 

6.78 

20 

25.7 

7.78 

1 

0.25 

102000 

0.017 

(1) 

2-L-4 

2.57 

0.13 

6.78 

20 

25.7 

7.78 

1 

1.0 

106000  ' 

0.018 

(1) 

2-L-4 

2.57 

0.13 

6.78 

20 

25.7 

7.78 

1 

4.0 

138000 

0.015 

(1) 

2-L-4 

2.57 

0.13 

6.78 

20 

25.7 

7.78 

1 

24 

178000 

0.019 

(1) 

5-L-ll 

2.0 

i 0.5 

5 

i 20 

20 

1 

10 

5 

1 

1 0 

24000 

0.014 

5-L-l! 

2.0 

0.5 

5 

20 

20 

10 

5 

1 

, 0.25 

i 48000 

i 0.017 

5-L-ll 

2.0 

0.5 

5 

20 

1 20 

10 

5 

1 

1.0 

, 42000 

0.021 

5-L-ll 

2.0 

0.5 

5 

20 

20 

10 

5 

1 ’ 

4.0 

52000 

0.018 

5-L-ll 

2.0 

0.5 

5 

20 

20 

10 

5 

1 

24 

68000 

0.022 

3-L-I6 

2.0 

0.5 

5 

4 

20 

10 

5 

5 

0 

48000 

0.016 

3-L-I6 

2.0 

0.5 

5 

4 

20 

10 

5 

5 

0.25 

87000 

0.010 

3-L-I6 

2.0 

0.5 

5 

4 

20 

10 

5 

5 

! .0 

103000 

0.017 

3-L-I6 

2.0 

0.5 

5 

4 

20 

10 

5 

5 

4.0 

139000 

0.029 

5-L-ll 

2.0 

0.5 

5 

4 

20 

10 

5 

5 

24 

568000 

0.012 

(1) 

(1)  Data  for  one  crack  tip. 
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TABLE  XVI 


RETARDATION  PARAMETER  DATA  FOR  2219-T85I  ALUMINUM  ALLOY  PLATE  IN  ROOM  TEMPERATURE 
DESICCATED  AIR  FOR  TENSION-COMPRESSION  LOAD  CLASS 


HOLD  TIME  EFFECTS 


Specimen 

No. 

S 

R 

U 

c 

s 

Time  ^ ^5' 
Hours 

a*,  Inch 

4-L-l 

2.0 

0.1 

9 

EH 

20 

10 

-25 

0 

8500 

0.029 

4-L-l 

2.0 

0.1 

9 

BkI 

20 

10 

-25 

0.25 

6500 

0.024 

4-L-l 

2.0 

0.1 

9 

-0.8 

20 

10 

-25 

1.0 

7500 

0.029 

4-L-l 

2.0 

0.1 

9 

-0.8 

20 

10 

-25 

24 

8500 

0.036 

4-L-I4 

2.5 

0.1 

9 

-1.0 

25 

10 

-25 

0 

26000 

0.051  (1) 

5-L-ll 

2.5 

0.1 

9 

-1.0 

25 

10 

1 

-25 

0.25 

14000 

0.041 

5-L-ll 

2.5 

0.1 

9 

-1.0 

25 

10 

1 

-25 

1.0 

14000 

0.040 

5-L-ll 

2.5 

0.1 

9 

-1.0 

25 

10 

1 

-25 

24 

12000 

0.038 

2-L-6 

3.21 

0.13 

6.78 

-1.0 

25 

7.78 

1 

-25 

0 

52000 

0.024 

2-L-6 

3.21 

0.13 

6.78 

-1.0 

25 

7.78 

1 

-25 

0.25 

45000 

0.042 

2-L-6 

3.21 

0.13 

6.78 

-1.0 

25 

7.78 

1 

-25 

1.0 

40000 

0.042 

2-L-6 

3.21 

0.13 

6.78 

-1.0 

25 

7.78 

1 

-25 

24 

54000 

0.057 

l-L-19 

2.5 

0.5 

5 

-1.0 

25 

10 

5 

-25 

0 

24000  - 

0.015  (1) 

4-L-l 

2.5 

0.5 

5 

-1.0 

25 

10 

5 

-25 

0.25 

24000 

0.016 

4-L-l 

2.5 

0.5 

5 

-1.0 

25 

10 

5 

-25 

1.0 

24000 

0.016 

4-L-l 

2.5 

1 0.5 

5 

-1.0 

25 

10 

5 

-25 

24 

24000 

0.026 

^-L-I3 

2.5 

0.1 

9 

-2.0 

25 

10 

1 

-12.5 

0 

29000 

0.056  (1) 

l-L-9 

2.5 

0.1 

9 

-2.0 

25 

10 

1 

-12.5 

0.25 

24000 

0.050 

l-L-9 

2.5 

0.1 

9 

-2.0 

25 

10 

1 

-12.5 

1.0 

22000 

0.053 

l-L-9 

2.5 

0.1 

9 

-2.0 

25 

10 

1 

-12.5 
I 

24 

20000 

0.060 

(I)  Data  from  basic  program. 
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TABLE  XVII 


PREDICTED  DELAY  CYCLES  USING  DIFFERENT  NUMERICAL  INTEGRATION  STEPS 
FOR  K|=I5,  K2=I0,  AND  KSI  /InT 


NUMBER  OF 
INTEGRATION  STEPS 
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G.  CLOSURE 

25 

7405 

2391 

4895 

50 

7411 

2391 

4895 

100 

7413 

2391 

4896 

250 

7413 
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4896 
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PLATES  SUBJECTED  TO  TENSION-ZERO  OVERLOADS  AND 

K,  = 10  ksi/iT^ 


117 


TABLE  XXI 

MEASURED  AND  PREDICTED  DELAY  CYCLES  USING  VARIOUS 


118 


TABLE  XXII 

MEASURED  AND  PREDICTED  DELAY  CYCLES  USING  VARIOUS 
RETARDATION  MODELS  FOR  2219-T851  ALUMINUM  ALLOY 
PLATES  SUBJECTED  TO  TENSION-COMPRESSION  OVERLOADS 
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SPECIMEN 

K 

K 
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NO. 

S 

R 

^1 

'^5 

G.  WHEELER 

G.  WILLENBORG 
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0.1 

15 

-15 

4500 

1944 

3653 
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1.5 

O.I 
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0.3 
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APPENDIX  A 

RETARDATION  MODEL  FORMULATIONS 


A-1 


APPENDIX  A 


This  appendix  contains  definition  and  discussion  of  the  three  retardation  models 
used  for  predicting  delay  cycles  for  the  cases  evaluated  experimentally.  In  order  of 
discussion,  the  models  were:  1)  generalized  Wheeler  model,  2)  generalized  Willenborg 
model,  and  3)  Grumman  closure  model. 


1 . GENERALIZED  WHEELER  MODEL 

The  first  retardation  model  was  proposed  by  Wheeler  (Al)  who  suggested  the  use 

of  a reduction  factor,  C , on  the  steady-state  (constant-amplitude)  crack-growth  rate  to 

P 

obtain  the  retarded  crack-growth  rate  under  spectrum  loading  as 


rda 

^dN 

sp 


/da 

ss 


and  C is  defined  as 
P 

m 
z 

(“|)  for  z < z* 

z*  c 

1 for  z ^2* 

c 


(Al) 


(A2) 


In  Equation  (A2),  the  quantity  z , equivalent  to  R in  the  original  Wheeler  paper, 

c y 

is  the  extent  of  the  current  load  interaction  zone  caused  by  the  current  maximum  applied 

cyclic  stress  intensity  factor,  K . The  quantity  z*,  which  is  equivalent  to  a - a 
' ' max  ^ p 

in  Reference  (Al),  is  the  difference  between  the  load  interaction  zone  due  to  a previous 
overload,  the  current  crack  growth  increment  since  that  overload,  Aa,  or 

z*  = Zqi^  - Aa  (A3) 

z*  is  also  the  load  interaction  zone  that  would  be  necessary  to  have  no  retardation. 

Wheeler  assumed  that  the  load  interaction  zone,  z,  was  equal  to  the  plastic 
zone  radius  created  under  plain-strain  condition,  i.e. 


A-2 


2 


1 X 

_ I ! max^ 
ys 


where  is  the  yield  strength  of  the  material. 


(A4) 


For  any  given  load  interaction  zone,  z*,  there  always  exists  a such  that 

Equation  (A4)  is  satisfied.  By  substituting  Equation  (A4)  into  Equation  (A3)  and  rearrang- 
ing the  results,  one  obtains 

1/2 


K*  = [1  - — ] 

max  max  ^OL 

where  the  maximum  stress  intensity  factor  generated  by  the  overload. 


(A5) 


The  steady-state  and  spectrum  crack-growth  rates  can  be  written  in  the  following 
general  forms: 


and 


& =c[f(AK,  R)l" 
ss 


sp 


(A6) 


(AT) 


respectively.  In  the  above  equations,  c and  n are  constants  for  a given  stress  ratio, 

R,  f(&K,R)  is  the  controlling  stress  intensity  function,  and  is  an  effective  stress 

intensity  range  under  variable-amplitude  loading.  The  function  f(AK,R)  has  different 
forms,  depending  upon  the  growth  rate  equation  used,  e.g., 

f(AK,R)=AK  (A8) 

if  Paris'  Equation  (A2)  is  used; 

m. 


f(AK,R)  = (K  -Kj- 
max  th  n 


H 


(A9) 


.AK 


if  the  Hall  Equation  (A3)  is  used; 

f(AK,R)  = (1  +0.6R)aK 
if  the  Grumman  Equation  (A4)  is  used;  and 


(AlO) 


A-3 


f(AK,R)  = K (1-RrW 


max 


(All) 


if  Walker's  Equation  (A5)  is  used, 


Substitution  of  Equations  (A2)  through  (AT)  into  Equation  (Al)  and  after  some 
simplification  lead  to 

-i2m/n 


AK 


eff 


K 


max 


max  ^OL 


172 


f(AK,R) 


= f(AK,R) 


ifK*  >K 
max  max 


if  K*  5 K 
max  max 


(A12) 


Experiments  involving  single  overloads  in  2024-T3  aluminum  and  Ti-6Ai-4V  (A6) 

indicate  that  there  is  a particular  value  of  overload  ratio,  say  S , such  that  when 
OL 

K /K  g crack  arrest  occurs*  This  particular  value,  S , is  called  the 

max  max  so  so 

"overload  shut-off  ratio"  and  its  corresponding  overload  is  called  the  "shut-off 
overload. " 


By  assuming  that  the  shut-off  overload  level  develops  the  local  stress  state  such 
that  the  effective  stress  intensity  range  just  equals  the  threshold  stress  intensity  range, 
AKth#  the  point  below  which  no  measurable  fatigue  crack  growth  occurs,  then,  from 
Equation  (A12),  one  obtains 


2m/ n 

-)  . f(AK,R) 


(A13) 


so 


from  which  yields  the  exponent  m as 


■ j AK^^ 
|f(AK,R)! 


m 


log  (5) 


SO  -J 


or 


n 

m - 2 


f(AK,R) 


log  aK 


th 


log  S 


so 


(A14) 


It  is  clear  that  the  exponent  m depends  not  only  on  the  specific  material 
being  used,  but  also  on  the  loading  subsequent  to  the  overload. 


A-4 


For  any  load  cycle  subsequent  to  the  overload,  having  crack  growth  increment 
since  the  overload,  ia,  one  can  calculate  its  effective  stress  intensity  range, 
from  Equations  (A12)  and  (A14),  and  its  corresponding  crack  growth  rate  from 
Equation  (AT). 


2.  GENERALIZED  WILLENBORG  MODEL 

Since  the  rate  of  propagation  of  a fatigue  crack  is  controlled  by  the  stress 
intensity  factor  at  the  crack  tip,  the  magnitude  of  the  stress  intensity  factor  is  a good 
indicative  of  the  extent  of  crack  tip  deformation.  Willenborg,  Engle,  and  Wood  (AT) 
then  proposed  a model  using  an  "effective  stress"  concept  to  reduce  the  applied  stress, 
and  hence  the  crack  tip  stress  intensity  factor.  To  help  identify  the  physical  implica- 
tions, Gallagher  (A8)  reformulated  the  Willenborg  et  al.  model  by  using  a residual  stress 
intensity  factor  concept  which  suggests  that  the  local  or  effective  stress  intensity  factor, 
be  calculated  from  the  equation 

(AI5) 


where  K is  the  stress  intensity  factor  associated  with  far-field  applied  load,  and 

Kq  is  the  residual  stress  intensity  factor  which  has  the  following  general  form: 

K 


K_  = , K^*-,  N„, , K , R,  a , N,  . . .) 

R max  min  OL  max  ys 


(A16) 


The  undefined  parameters  in  the  above  equation  are 

= minimum  K generated  by  the  overload 

mm 


Nql  ^ number  of  overloads  applied 


N 


number  of  cycles  of  K 


max 


The  equivalent  residual  stress  intensity  factor  for  the  Willenborg  et  al.  model  is 

1/2 


(1  - -^) 

R max 


K = K*  - K 
max  max  max 


(A17) 


where  's  the  overload  created  load  interaction  zone,  which  is  approximately  equal 
to  the  overload  plane  stress  plastic  zone  radius 
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2 


OL  y 


,OL 


max 


OL 


2tr 


/s 


and  Aci  !s  an  increment  of  crack  growth  into  the  overload  load  Interaction  zone  since 
the  application  of  overload.  Following  an  overload,  Aa  is  approximately  zero  and  the 
residual  stress  intensity  factor  is  maximum. 


Equation  (AIT)  has  been  found  in  error,  since  crack  growths  have  been  observed 
after  single  overload  application  where  the  overload  ratio,  which  is  defined  as 


K 


S = 


OL 

max 


K 


max 


is  greater  than  or  equal  to  2.0,  the  condition  for  which  the  Willenborg  et  al.  model 

predicts  zero  crack  tip  effective  stress  intensity  factors.  - Gallagher  and  Hughes  (A9) 

suggested  that  the  actual  residual  stress  intensity  factor  may  be  proportional  to  the 

W ^ 

Willenborg  et  al.  residual  stress  intensity  factor,  , and  introduced  the  factor  0 
such  that 

W 

K|^  -0K|^ 

Using  Equations  (A17)  and  (A20),  the  effective  maximum  and  minimum  stress 
intensity  factors  can  be  written  as 

1/2 


and 


(K  ) K "0 
max  --  max  ^ 
eff 


(K  . ) = K . "0 

mm  -r  mm  ^ 
ett 


K°'-  (1  - 

max  ^OL 


(1  - ) 
max  ^OL 


K 


max 


1/2 


- K 


max 


or  (K  . ) , whichever  is  larger. 


mm 


eff 


cf 


where  (K  . ) is  the  cut-off  value  of  the  minimum  effective  stress  intensity  factor. 


min 


eff 


Normally,  ) Is  set  equal  to  zero.  When  the  overload  ratio  is  greater  than  3.0, 


eff  , 
cf 

it  has  been  found  that  (K  . ) = -0,3  (K  ) will  give  a better  prediction. 

m in  max 

eff  r eff 

cf 


(A18) 


(A19) 


(A20) 


(A21) 


(A22) 


The  effective  stress  intensity  range  and  effective  stress  ratio  are  then  calculated 
from  the  following  equations: 

(aK)„  = (K^  ) “(K  ) (A23) 

ett  max  min 

ett  ett 


and 


(K  :_) 


min 
max 


eff 


(A24) 


eff 


By  assuming  that  the  shut-off  overload  level  develops  a local  (effective)  stress 

intensity  factor  condition  such  that  no  growth  is  induced,  Gallagher  et  al.  (A9)  set  the 

maximum  local  stress  intensity  factor,  , equal  to  the  maximum  fatigue  threshold 

stress  intensity  factor,  since  approximately  constant  for  negative 

stress  ratios,  (A10)  Immediately  following  the  shut-off  overload  which  produces  no  growth, 

Aa  = 0 and  (K^^qx^  ~ ^^max)  ' solve  Equation  (A21)  for  0 subjected  to  these 

^ eff  TH 

assumptions  as 


0= 


1 - (K  ) /K 

max  -L,  max 
I n 

___ 

so 


(A25) 


For  any  load  cycle  subsequent  to  the  overload,  one  can  calculate  its  effective 
stress  intensity  range,  from  Equations  (A21)  through  (A25),  and  its  corresponding 

crack  growth  rate  from  any  crack  growth  rate  equation  by  setting  aK  = ^ ~ ^eff  * 


3.  GENERALIZED  CLOSURE  MODEL 

In  1968,  Elber  (All)  observed  that  the  surface  of  the  fatigue  crack  closed  during 

fatigue  crack  propagation  tests  under  cyclic  tension  loading  (tension-tension  load  cycle) 

on  aluminum  alloy  sheet  specimens.  Based  on  results  of  tests,  he  suggested  that  the  tip 

of  a fatigue  crack  closes  while  specimens  were  still  subjected  to  a tension  loading,  and 

that  the  fatigue  crack  propagates  only  during  that  portion  of  the  load  cycle  in  which  it 

was  fully  open.  Assuming  that  the  crack  closure  stress  is  not  different  from  the  crack  open 

stress,  then  according  to  Elber*s  hypothesis,  the  effective  stress  range,  is  equal 

to  the  difference  between  maximum  stress,  a , and  the  closure  stress,  a , i.e., 

max  c 
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(A26) 


or 


A cr  r r = cr  - a 

err  max  c 


ACT  rr  - CT  (1  - C J 
erf  max  f 


where  is  the  closure  factor  defined  as 


(A27) 


max 


The  closure  stress,  in  general  is  a function  of  applied  stress  ratio,  previous 
applied  load  history  and  any  other  pertinent  parameters. 


Bell,  et  al,,  (A4)  uses  the  Elber  crack-closure  concept  and  modifies  the  Paris 
equation  as 


da 

dN 


- c 


max 

V 


o 


-v^ 


(A28) 


for  the  case  of  a through-crack  in  an  infinite  sheet.  The  coefficients  c and  n in 

Equation  (A28)  have  the  same  values  that  would  be  determined  by  using  the  Paris 

equation  to  fit  crack-growth  rate  data  plotted  against  the  actual  applied  stress  intensity 

range.  is  the  closure  factor  at  R = 0. 

o 

Equation  (A28)  can  be  modified  for  any  other  geometric  configuration  by  intro- 
ducing the  total  correction  factor  to  account  for  all  the  appropriate  stress  intensity 
magnification  factors,  as 


da 

or 

- 0 



max 

c 

dN 

c 

(1  - 

Cf) 

o 

(A29) 


Equation  (A29)  can  also  be  expressed  in  terms  of  stress  intensity  factors  as 
da 


dN 


c 


K - K 
max  c 


1 - C 


f 


(A30) 


where  K is  the  closure  stress  intensity  factor, 
c 


Following  a similar  approach,  the  Hall  and  Walker  crack  growth  rate  equations 
may  be  written  in  the  closure  form  as 

da  _ 

Tn 

and 

da  ^ 

dN 

respectively. 

a.  Closure  Model  for  Constant  Amplitude 

The  constant-amplitude  crack  growth  rate  equation  for  2219-T851  aluminum 
was  given  as  (A4) 

-^  = 1.96  X 10“^  [(1  +0.6R)  forO  s R s 0,5  (A33) 

-9 

by  letting  c = 1 .96  x 10  and  n = 3.34  in  Equation  (A30),  and  solving  Equations  (A30) 
and  (A33)  for  the  closure  factor  leads  to 

K 

= ^ = 1 - (1  - Cf  )(1  + 0.6R)(1  - R)  (A34) 

max  o 

However,  based  on  the  test  data  generated  under  various  constant  stress  ratio 
(-1  S R § 0.7),  the  closure  factor  was  later  chosen  as  (A4) 

Cj,  = Cj,  + (Cj.  - Cj.  )(1  + R)^  (A35) 

-1  o -1 

For2219-T851  aluminum,  C^.  =0.347,  C^.  =0.40,  and  P = 3.93. 

-1  o 

Equation  (A35)  was  also  proposed  to  account  for  the  effects  of  negative  stress 
ratios  and  compressive  stresses. 

For  constant-amplitude  load  cycles,  Equation  (A35)  is  used  to  calculate  the  closure 
factor  Cj..  The  corresponding  closure  stress  and  its  crack  growth  rate  are  computed  from 


m. 


(K  - K . ) 

max  th 


K - K 
' max  c 

O J 


H 


(A31) 


- n. 


K - K 
max  c 

1 - C. 


w 


(A32) 


A-9 


Equations  (A27)  and  (A30),  respectively,  where  in  Equation  (A30)  c = 1.96  x 10  , 

n = 3.34,  and  = 0.40. 
o 

b.  Closure  Model  for  High- Low  Loading  Sequence  Having 
Constant  Minimum 

Let  and  be  the  constant-amplitude  closure  stresses  of  high  stress 
cycles,  Oy  and  lower  stress  cycles,  respectively.  Then,  fora  high-low  loading 
sequence  having  a constant  minimum  stress  level,  the  crack  growth  during  the  lower 
loading  cycles  will  be  retarded  before  it  changes  to  the  steady-state  growth.  During 
the  retarded  period,  the  closure  model  assumes  that  the  closure  stress,  a^,  varied  accord- 
ing to  the  following  equation: 

. B 

a = a - (a  - a ) (—)  for  0 5 Aa  ^ z (A36) 

c c^  ' C2  z 

where  Aa  is  the  crack  growth  since  the  stress  change,  z is  the  load  interaction  zone 
caused  by  Oj , and  B is  an  empirical  exponent.  A process  of  trial  and  error  was  used 
to  obtain  the  values  of  the  parameters  for  2219-T851  aluminum  as 

=0.4 

o 

z ^ plane  stress  plastic  zone  radius 
B ^1.0 


c.  Effect  of  Minimum  Stress  on  Closure 


If  the  minimum  stress  of  applied  load  cycles  decreases  from  one  load  level 

to  the  second  load  level,  regardless  of  the  maximum  stresses,  the  closure  model  assumes 

that  step-wise  change  in  closure  stress  level  occurs  during  the  first  load  cycle  after  the 

load  change.  The  closure  factor  immediately  after  the  application  of  the  new  lower 

minimum  stress,  a . , is  calculated  as  C*  from  Equation  (A35)  using  R = a • /a, 

mm2  f \ mm2 

and  its  corresponding  closure  stress  is  calculated  from  Equation  (A27)  as  a*  = a,  C'r* 

C ] IT 

Starting  from  this  value,  the  closure  stress  varies  as  described  by  Equation  (A36).  The 

same  concept  can  be  extended  to  the  case  where  the  minimum  stress  is  negative,  i.e. 

0*.^*  <0,  For  the  case  where  R < -1,  the  current  closure  model  treats  it  the  same  as 

mm2 

R = -1. 
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d.  Increasing  and  Decreasing  Closure  Level 


To  use  the  closure  model,  it  is  necessary  to  define  whether  the  closure  ^ 

level  is  increasing  or  decreasing  during  variable-amplitude  loading.  The  current  closure 
stress,  which  is  defined  in  the  value  which  existed  prior  to  a load  change,  may  be 
modified  by  overload  or  minimum  stress  consideration  immediately  after  the  application 
of  the  new  loading.  If  the  current  closure  stress  is  greater  than  the  expected  closure 
stress,  which  is  defined  as  the  stabilized  value  associated  with  the  new  loading,  a 
decreassing  closure  situation  exists,  and  the  behavior  of  the  closure  stress  is  defined  by 
Equation  (A36  ). 


Conversely,  when  the  closure  level  is  increasing,  i.e.  the  current  closure  stress 
is  less  than  the  expected  closure  stress,  say  a , then  the  closure  stress  is  calculated 
from  the  following  equation  for  2219-T851  aluminum: 


(A37) 


or  1 whichever  is  smaller 

where  is  the  number  of  consecutive  overload  cycles  after  a load  change.  Accord- 

ing to  Equation  (A37 ),  13  cycles  of  the  new  loading  are  required  to  regain  a stabilized 
closure  value  for  2219-T851  aluminum. 
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